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1 .  INTRODUCTION 


Large  amounts  of  global  cloud  cover  data  in  the  form  of 
visible  and  infrared  satellite  imagery  are  gathered  daily  by 
ground  stations  throughout  the  world.  These  data  are  readily 
accessible  and  available  for  exploitation.  Nevertheless,  the 
extraction  of  useful  information  (cloud  type,  fractional  cov¬ 
er,  precipitation  amounts)  from  satellite  data  is  not  keeping 
pace  with  the  rate  and  volume  at  which  these  data  are  collect¬ 
ed  . 

Although  valuable  when  practiced  by  a  skilled  analyst, 
current  subjective  satellite  imagery  cloud  classification 
techniques  are  time-consuming  and  costly.  Hard  copies  are 
expensive  to  produce  and  distribute  on  a  large-scale  basis. 
Accurate  human  analysis  of  all  available  data  is  not  feasible, 
and  such  analyses  may  lack  consistency  of  interpretation  among 
photoanalysts . 

The  abundance  of  satellite  data  collected  by  Air  Force 
Global  Weather  Central  (AFGWC)  has  sparked  interest  in  the 
development  of  efficient  real-time  imagery  processing  me¬ 
thods  which  provide  meteorological  information  quickly  and 
accurately.  AFGWC  currently  utilizes  automated  cloud  analysis 
and  classification  procedures  which  determine  3-dimensional 
global  cloud  cover  distributions  using  a  25  nautical  mile  (nm) 
grid  spacing.  This  process,  called  3D  Nephanalysis  (3DNEPH), 
uses  imagery  data  from  the  Defense  Meteorological  Satellite 
Program  (DMSP)  along  with  surface,  radiosonde,  and  aircraft 
data.  DMSP  satellites  are  equipped  with  very  high  resolution 
or  fine  mode  (1/3  nm)  visible  and  infrared  sensors.  Nonethe¬ 
less,  satellite  tape  recorders  and  AFGWC' s  computer  limita¬ 
tions  restrict  3DNEPH  satellite  input  to  smoothed  (3  nm) 
imagery  for  global  applications.  A  decrease  in  the  quantity 
of  input  inevitably  decreases  the  quality  of  output,  especi¬ 
ally  for  areas  where  conventional  data  are  sparse  or  not 
available . 

Hence  there  is  a  need  for  improving  real-time  satellite 
imagery  processing  capabilities.  The  incorporation  of  fine 
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mode  satellite  imagery  is  fundamental  for  fulfilling  the 
need,  even  though  significantly  larger  amounts  of  data  will 
have  to  be  processed.  Past  and  reliable  analytical  proces¬ 
sing  techniques  are  reasonable  and  wise  approaches  to  util¬ 
izing  fine  mode  imagery  data.  Techniques  must  be  able  to 
reduce  large  amounts  of  satellite  imagery  data  while  retain¬ 
ing  as  much  information  as  possible  about  any  clouds  in  the 
imagery . 

Discrete  Fourier  transform  power  spectrum  analysis  is 
the  basis  of  one  such  technique  which  has  been  studied  and 
is  described  in  this  report.  An  algorithm  was  developed  as 
an  alternative  to  that  part  of  3DNEPH  which  determines 
cloud  types.  For  this  study  the  algorithm  was  updated  and 
modified  for  use  with  2/3  nm  DMSP  visible  and  infrared  data. 
It  classifies  a  25x25  nm  imagery  sample  as  one  of  nine  pos¬ 
sible  cloud  types  in  less  than  one  second  on  the  Air  Force 
Geophysics  Laboratory  (AFGL)  CDC  6600  computer. 

This  report  presents  a  detailed  mathematical,  analyti¬ 
cal,  and  computational  description  of  this  cloud  classifica¬ 
tion  technique.  It  also  discusses  the  significance  of  its 
results,  and  offers  some  ideas  for  further  refining  and  im¬ 
provement  of  the  algorithm. 

2 .  DATA 

The  spectral  analysis  routine  uses  visible  and  infra¬ 
red  data  from  a  DMSP  satellite  in  a  450  nm  sun- synchronous 
orbit.  In  particular,  imagery  from  several  Mediterranean- 
Middle  East  orbits  was  studied.  Figures  la  and  lb  show 
samples  of  visible  and  infrared  imagery  from  orbit  7399+9 
of  DMSP  vehicle  F-l .  The  Sinai  peninsula  (A)  is  located  in 
the  lower  center  of  each  image.  To  the  west  of  the  Red  Sea 
and  the  Gulf  of  Suez,  the  Nile  River  can  be  seen. 

Each  data  value,  or  pixel,  represents  the  brightness  or 
temperature  of  a  2/3  nm  region  of  the  earth.  2/3  nm  DMSP 
data  are  unique  to  this  study.  DMSP  satellites  are  equip¬ 
ped  with  1/3  nm  sensors.  At  AFGWC ,  the  brightness  value  at 
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Figure  la.  Visible  imagery  from  DMSP  orbit  #7399+9  (February  1978),  showing  only  every 
fifth  picture  element  of  every  fifth  line  of  2/3  nm  imagery.  A  marks  Sinai  Peninsula. 
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each  point  of  the  2/3  nm  imagery  was  obtained  by  averaging 
four  pixels  taken  from  the  original  1/3  nm  data. 

The  size  of  a  sample  that  is  classified  by  this  algor¬ 
ithm  is  37x37  pixels.  It  represents  a  25x25  nm  region  of  the 
earth,  which  is  3DNEPH's  horizontal  grid  spacing,  since  each 
pixel  represents  2/3  nm. 

The  Man-computer  Interactive  Data  Access  System  (Mc- 
IDAS)  at  AFGL  was  used  to  display  the  DMSP  imagery  and 
store  its  digital  values.  McIDAS  is  capable  of  storing  24 
500  x  672-pixel  arrays  per  user  terminal.  Thus  the  data  for 
up  to  12  visible  and  12  infrared  images,  each  one  covering  a 
333  x  448  nm  area  of  the  earth,  can  be  readily  accessed. 

(See  Figures  2a  and  2b.)  The  invoking  of  a  McIDAS  save  com¬ 
mand  enables  a  user  to  search  interactively  for  interesting 
cloud  samples  within  a  displayed  DMSP  image.  Once  such  a 
sample  is  found,  the  user  subjectively  evaluates  the  cloud 
type  in  the  image  sample,  along  with  its  background,  for  veri¬ 
fication.  This  information,  in  addition  to  the  sample's 
visible  and  infrared  satellite  data,  is  then  written  on  a 
sample  save  tape.  This  process  is  continued  until  the  user 
feels  an  adequate  cloud  imagery  sample  set  for  use  by  the 
automated  classifier  has  been  collected. 


Figures  2a  and  2b  show  the  visible  and  infrared  DMSP 
imagery  of  an  area  of  northern  Egypt  and  the  northwestern 
region  of  the  Sinai.  The  Nile  delta  is  clearly  visible  in 
the  center  of  Figure  2a.  The  darker  area  in  the  lower  cen¬ 
ter  of  the  image  is  El  Faiyum,  a  topographical  depression 
irrigated  by  the  waters  of  the  Nile.  The  Gulf  of  Suez  is 
in  the  lower  right.  A  line  of  cumuliform  clouds  lies  to 
the  west  of  the  Nile  (left-center  of  the  image)  covered  by 
a  thin  cirrostratus  shield,  hardly  detectable  in  the  visible 
but  highly  detectable  in  the  infrared.  A  25  x  25  1®  cloud 
imagery  sample  used  by  the  spectral  analysis  routine  covers 


the  size  of  the  area  covered  by  Figures  2a  and  2b. 
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Figure  2a.  Visible  imagery  from  DMSP  orbit  #7399+9  (February  1 978).  This  image  is 
original  2/3  nm  resolution.  Area  shown  is  the  Nile  delta. 


3.  DESCRIPTION  OF  THE  CLASSIFIER 

The  spectral  classifier  technique  consists  of  two  main 
parts i  the  discrete  Fourier  transform  and  a  Bayes  decision 
rule.  The  discrete  Fourier  transform  of  a  cloud  imagery  sam¬ 
ple  is  first  computed,  since  the  transform  contains  informa¬ 
tion  about  the  shapes  and  sizes  of  clouds  present  in  the 
imagery  field.  The  transform  is  then  reduced  in  size  to  a 
power  spectrum,  which  is  a  vector  whose  components  are  nor¬ 
malized  averages  of  the  coefficients  found  in  the  transform. 
Then  a  Bayes  discriminant  function  decides  the  cloud  type  of 
the  satellite  imagery  sample  using  the  information  contained 
in  the  power  spectrum.  The  discriminant  function  is  derived 
from  the  multivariate  normal  density  function. 

In  the  remainder  of  the  section,  a  mathematical  formali¬ 
zation  of  this  technique  is  developed  in  order  to  provide  a 
better  understanding  of  the  spectral  classifier.  Symbols  are 
defined  in  the  Glossary  of  Terms,  page  74. 

The  continuous  Fourier  transform  has  been  a  standard 
analytical  tool  for  solving  many  types  of  pattern  recognition 
problems.  However,  until  recently,  applications  of  the  dis¬ 
crete  Fourier  transform  (DFT)  have  been  limited  because  of 
the  large  amounts  of  computational  time  required  to  calculate 
it.  With  the  development  of  the  fast  Fourier  transform,  an 
algorithm  which  computes  the  DFT  quickly  and  efficiently, 
the  DFT's  use  in  problem  solving  has  greatly  increased. 

3.1  The  Discrete  Fourier  Transform 

A  real-valued,  continuous  function  F(x,y)  of  two 
variables  defined  on  the  region  where  0<x<M,  0<y<N, 

M  and  N  positive  integers,  can  be  expressed  as  a  unique 
linear  combination  of  products  of  sines  and  cosines « 
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for  (x,y)C^?,  and  where  A  ,  B  ,  C  ,  and  D  C  R. 

m,n  m,n  m,n  m,n 

This  is  called  the  double  Fourier  series  expansion  of  F 
on  m.  This  representation  extends  F  periodically  for  all 
(x,y)€R  .  An  alternate,  shorter  form  is: 

F  (x  ,y)  =  ^  ip  2  -n\£i  +  nJ)  iji  CC 

^ ^m,n  ’  m,n  ( 

m=-°°  n=-°° 


where  M  N 

^m.n  =  MN  y^^x,y^e"2  +  ~”^dydx»  Vm.nCZ.  (3) 

Note  that  i|i0  Q  is  the  average  value  of  F  over  a.  Equation 
(3)  is  known  as  the  double  Fourier  transform  of  F  on  M. 

There  are  infinitely  many  harmonics  in  F,  as  can  be 
seen  by  looking  at  the  upper  and  lower  bounds  of  the  ^  /s 
in  (1)  and  (2).  Infinitely  many  coefficients  ♦  can  be 
solved  for  by  (3)  since  F  is  continuous  on  tfl.  Equations 
(1)  and  (2)  are  exact  representations  of  F(x,y)  on  @1. 

Let  ^  =  { (  i  ,  j  ) |  (  i  ,  j  ) €<s^?and  i,jCZ},  where  Z  is  the  set 
of  integers.  Consider  the  discrete  function 

F(i,j)  =  F  (  i  ,  j  )  V  (  i ,  j  ) 

A 

F  can  be  extended  to  the  entire  lattice  of  points  (i,j)  on 
2 

the  plane  R  such  that 

F  (  i+kM,  j+lN)  =  F(i,j)  for  any  k.lCZ, 

A  A 

that  is,  F is  periodic  on  Z x  Z .  F  can  be  expressed  as  a 
finite  series  in  the  form: 

[f]  [I] 

*(1> ^  ^  j  Vncos  ^  cos  *  ®m,ncos 
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where  A  .  B  ,  C  D  „CR,  This  is  called  the  discrete 
th  i  n  in  i  n  m  •  n  m  i  n  ^ 

double  Fourier  series  expansion  of  F  on  kK  .  Its  shorter, 
complex  form  is: 


F(i,j)  = 


%  „ 


cc>  (s) 


where  [x]  is  the  greatest  integer  less  than  or  equal  to  x, 

x  positive,  or  the  smallest  integer  greater  than  or  equal  to 

x,  x  negative.*  The  $  are  given  by  equation  (6): 

M  N  m,n 

*m,n  -  ®  'C1<SM  *  m.nCZ 


such  that 


Iff] 


<  m  <  l^l .  and 


<  n  <  . 


(6) 


Note  that  <ji0  Q  is  the  average  value  of  F  over^. 
Equation  (6)  is  known  as  the  discrete  double  Fourier  trans- 

^  A  1  —  ■  —  -  ■ 

form  of  F  on 

The  coefficients  in  any  Fourier  series  expansion  of  a 

/“V 

discrete  function  F  are  the  amplitudes  of  the  trigonometric 
.  o 

functions  they  precede.  '  Each  coefficient  in  this  expansion 
contains  information  pertaining  to  the  size  of  a  discrete  har- 
monic  that  is  present  in  F. 

/\ 

Since  it  is  discrete,  F  contains  only  a  finite  number 
of  harmonics.  There  is  a  limit  to  the  number  of  unique  co¬ 
efficients  £  that  can  be  determined.  This  statement  is 

in  1  n 

proved  in  Appendix  B. 

In  a  discrete  sense,  the  smallest  resolvable  nontrivial 
harmonic  (wave)  in  either  the  x-  or  y-direction  needs  at 
least  two  points  to  establish  itself,  so  that  the  largest 

A 

wavenumber  F  can  possibly  contain  in  the  x-direction,  given 
M  points  in  the  x-direction,  is  jjfJ  ,  and  in  the  y-direction, 
given  N  points  in  the  y-direction,  is  j¥l  .  It  is  impossible 


*See  Appendix  A  for  the  relationship  between  (5)  and  (4). 

2 

This  is  easy  to  see  in  equation  (4).  In  equation  (5). 
remember  e  =  cos  0  +  isine  ,  i  *  /^1 . 
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arc  ~ore  than  \-x\  times  m  the  y-c:rection. 

Thus,  m  equation  (4),  there  are  ~  •  —  A  * s,  B  ’s, 
a  ^  MN  22m, n  m,n 

Cm  n*s,  and  n*s,  or  4  •  — ^  =  MN  distinct,  nontrivial 
solvable  Fourier  coefficients.  Likewise,  in  (5)»  there  are 
also  MN  distinct  i  n's‘  §eneral  only  MN  coefficients  can 

be  computed  when  F  is  represented  by  MN  values.-^ 


The  numbers  ^  and 


are  called  the 


wave- 


numbers  of  F  on  m  the  x-  and  m  the  y-directions , 
respectively. 


3.1.1  The  Significance  of  Fourier  Coefficients 


Each  coefficient  in  the  Fourier  expansion  of  a  discrete 
function  contains  information  pertaining  to  the  amplitude  of 
a  discrete  harmonic  that  is  present  in  that  function.  For 

.  .  A 

example,  if  there  exists  no  wave  m  F  with  wavenumber 
m  =  i,  j--|  <  i  ^  £|J  in  the  x-direction  and  n  =  j,  |~J  <  j 

<  in  the  y-direction,  then  the  Fourier  coefficients  in 
(4)  and  in  (5)  with  the  subscript  (wavenumber)  "i,j"  are 
equal  to  zero. 

The  harmonics  of  a  Fourier  expansion  are  generated  about 
the  mean  of  the  function  which  the  expansion  represents,  as  is 
shown  in  Figure  3* 

The  Fourier  series  representation  of  the  discrete  func¬ 
tion  of  one  variable  shown  in  Figure  3a  is  (set  j ,  n  =  0  in 
(4)): 

y  =  0.75  -  0.25  cos  f x, 

for  x  =  1,2,3,. ..,8.  Note  how  the  cosine  wave  (dashed  curve) 
of  wavenumber  4  and  amplitude  0.25  (A^  =  -0.25)  lies  about 
the  line  y  =  0.75.  which  is  the  mean  of  y. 

In  two  dimensions,  the  mean  of  a  function  is  not  a  line 


3  A  ,  ,  , 

^Just  as  f  is  periodic,  so  is  its  transform.  A  function 

A  ^ 

F(i,i)  of  space  is  transformed  into  a  function  ^  of 

m ,  n 

wavenumber . 


a.  b. 

Figure  3.  A  discrete  cosine  wave  y  of  amplitude  0.25  which  lies 
about  its  mean  O.75.  y  is  assumed  periodic  outside  of  the  region 
shown  here.  Note  how  the  "pattern"  established  in  the  first  two 
points  (x  =  1  and  x  =  2)  occurs  four  times  over  the  period  of  y, 
thus  making  A^  nonzero.  Figures  3a  and  3b  are  Fourier  transform 
pairs:  AQ  =  0.75.  and  A^  *  -0.25* 


Figure  4.  A  simple  discrete ,  ^finite  function  <£( i ,  j )  of  two  var¬ 
iables  defined  on  the  region  £0  ■  { ( i ,  j )  I  1  s  i  s  1  i 
and  j  integers}.  O(i.j)  is  assumed  periodic  outside  ofs#. 


but  a  plane  of  constant  distance  A^  ^  (which  is  equal  to 
from  the  x-y  plane.  The  harmonics  present  in  a  function  of 
two  variables  take  on  various  values  above  and  below  its 
mean  plane. 

Consider  the  discrete  function  G  of  two  variables  and 
its  Fourier  coefficients.  Let  M  =  N  =  4,  and  let  G  (i,  j)  be 

A 

defined  as  in  Figure  4a.  Computing  the  transform  of  G  for  the 

form  (4),  we  find  A1  .  =  3  and  all  other  15  coefficients  are 

1,1  £  . 

identically  zero.  This  is  true  because  G  is  a  discrete, 
pure  harmonic  of  wavenumber  m  =  1 ,  n  =  1 .  By  looking  at 
Figure  4b  it  is  more  easily  seen  that  there  is  in  fact  a 
discrete,  two-dimensional  harmonic  in  G  that  occurs 
exactly  once  horizontally  (m  =  1)  and  exactly  once  verti- 

A 

cally  (n  =  1).  Instead  of  defining  G  as  a  discrete  set  of 
sixteen  ordered  triples  (i,j,G(i,j)),  a  more  compact  ex¬ 
pression  in  the  form  (4)  would  be: 

,  .  /s  2ir  •  1  •  i  2tt  •  1  •  i 

G(i,j)  =  A1^1  cos  — £ — -  cos  — *L 

=  3  cos  ^  cos  ,  (?) 


where  i, j  C  Z  such  that  Is  i  s  4,  Is  j  s  4. 

A 

Now  suppose  we  define  a  new  discrete  function  H( i , j ) 

A 

as  in  Figure  5a.  Transforming  this  function  we  get  A„  , 

=  3,  and  all  other  15  coefficients  equal  zero.  His  a  pure 
"cosine-cosine"  wave  of  wavenumber  m  =  2 ,  n  =  2 .  As  can 
be  seen  in  Figure  5"b,  a  discrete  harmonic  occurs 
twice  in  the  horizontal  (m  *  2)  and  twice  in  the  vertical 
(n  =  2)  in  the  principal  part  of  its  domain,  and  has 
amplitude  3 ! 


H  (  i ,  j  ) 


cos 


2tt  -2 « i 
4 


cos 


;  2  •  j 
4 


=  3  COS  IT  i  COS  IT  j  .  (8) 

Equation  (8)  is  in  the  form  (4)}  i  and  j  are  as  defined 
for  G. 

As  was  seen  in  the  previous  two  examples,  the  Fourier 
transform  of  a  discrete  function  of  two  variables  is  sensi- 
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Figure  5-  A  simple  discrete,  finite  function  H(  i , j )  of  two 
variables  defined  on  the  region  ££>  as  defined  in  Figure  4. 
H(i,j)  is  assumed  periodic  outside  of  0. 
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tive  to  the  shapes  (wavenumbers)  and  sizes  (amplitudes)  of 
the  harmonics  in  that  function.  With  this  in  mind,  consider 
the  function  F  to  be  a  finite  set  of  visible  or  infrared 

A  I  ^ 

satellite  data  for  a  particular  region  on  the  earth.  F 
is  certainly  a  discrete  function  of  tvj£  variables.  As  long 
as  it  is  assumed  periodic  outside  of  ,*5.  its  transform  can 
be  computed.  The  coefficients  of  this  Fourier  transform  can 
provide  useful  information  pertaining  to  the  shapes  and  sizes 
of  any  clouds  present  in  F  over  M. 

This  information  is  the  basis  of  the  spectral  classi¬ 
fier  technique.  The  DFT  of  a  37  x  37-pixel  array  F  (25  x 
25  nm)  containing  the  satellite-sensed  visible  or  infrared 
brightnesses  of  a  cloud  sample  is  computed.  The  coefficients 
thus  obtained  may  then  be  examined  to  determine  the  texture 
and  patterns  of  the  clouds  in  the  imagery  sample. 

3.1.2  The  Fast  Fourier  Transform 

The  DFT  of  a  function  can  be  computed  from  equations 
like  (6),  but  for  large  M  and  large  N  this  is  a  slow  opera¬ 
tion.  It  is  then  more  rapidly  calculated  by  the  algorithm 
which  is  known  as  the  fast  Fourier  transform  (FFT). 

The  FFT  is  a  more  efficient,  quicker  way  of  computing 
the  coefficients  found  in  the  complex  Fourier  expansion  of 
a  discrete,  finite  function.  For  example,  the  FFT  algorithm 
calculates  the  Fourier  transform  of  a  one-dimensional  array 
of  1,024  values  over  200  times  faster  than  the  straight¬ 
forward  method.  Moreover,  the  output  of  an  FFT  is  inter¬ 
preted  in  exactly  the  same  manner  as  the  output  of  a  DFT. 

Several  restrictions  are  placed  on  any  data  entered 
into  the  FFT.  As  is  the  case  with  the  DFT,  the  number  of 
coefficients  in  the  transform  of  a  function  must  equal  the 
number  of  values  defining  that  function.  The  input  data 

^In  other  words,  a  finite  set  of platitude-longitude  grid  points. 
'’We  are  only  interested  in  f  on  and  nowhere  else,  so  this 
assumption  has  no  adverse  consequences. 
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must  represent  equally-spaced  values  of  a  function  in  the 
space  domain,  and  the  output  must  represent  equally-spaced 
coefficients  of  its  transform  in  the  wavenumber  domain.  The 
input  function  and  its  transform  are  assumed  to  be  one  com¬ 
plete  cycle  of  discrete  periodic  functions.  The  coefficients 
generated  by  the  FFT  are  not  normalized;  for  an  M  x  N  input 
array  they  appear  multiplied  by  M*N. 

In  general,  these  restrictions  are  minimal  and  the  FFT 
is  a  valuable  and  powerful  computational  tool  when  the  DFT 
of  a  large  function  of  two  variables  must  be  calculated. 


3.1.3  Annular  Integration 

/\ 

The  transform  of  a  37  x  37-pixel  array  F  contains 

37  ■  37  =  1.369  coefficients.  Analyzing  each  coefficient 

separately  would  be  a  painstaking,  time-consuming  approach. 

Therefore,  the  number  of  these  coefficients  must  somehow  be 

reduced,  while  at  the  same  time  retaining  as  much  information 

/\ 

about  the  clouds  in  the  imagery  sample  F  as  possible.  The 
data  reduction  technique  used  by  this  classifier  is  called 
annular  integration  and  is  described  here. 

The  square  root  of  the  sum  of  the  squares  of  the  real 
and  imaginary  parts  of  a  Fourier  coefficient  $  (6)  is 

defined  as  power ; 

P 

m,n 


Just  as  a  Fourier  coefficient  is  a  measure  of  how  far  a 
harmonic  deviates  from  the  mean  of  a  function  (amplitude), 
so  is  power,  since  it  is  defined  in  terms  of  the  components 


of  a  Fourier  coefficient.  The  set  of  powers  of  all  the  co- 

/V 

efficients  of  the  transform  of  a  function  F  is  called  the 


Fourier  spectrum  of  F. 


A  power  can  be  calculated  for  every  (m,n)  of  the  wave- 

number  plane  for  which  there  is  a  ^  .  These  terms  can  be 

m,n 

placed  onto  the  m-n  plane  by  mapping  each  Pm  to  the  point 
(m,n).  Then,  the  powers  can  be  sorted  into  annular  bands 


whose  components  are  the  same  p-distance  from  the  origin,  where 
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^  -v.  r*."\ 


do) 


p  -  [JJTT,  0.5],  for  ff]<m<|],  ana  [ff]<n<[f]. 
and  where  [  ]  is  as  previously  defined. 

Given  an  M  x  N  array  of  coefficients,  the  largest  number 
m  can  be  is  Jj[J  ,  and  the  largest  number  n  can  be  is  |j|J . 
Therefore  the  annular  band  with  the  largest  "radius"  P 


max 


(ID 


(12) 


must  be  (from  (10)): 

P  ! 

Kmax 

It  is  obvious  that 
pmin 

For  each  p,  every  power  term  in  annular  band  p  is 
squared,  and  these  squares  are  summed  together.  Finally, 
they  are  averaged,  this  average's  square  root  is  taken,  and 
then  it  is  normalized.  This  number  is  called  a  normalized 
average  amplitude  (NAA),  and  is  given  by: 


+  [f]2  ♦  0.5]. 

[/°2  ♦  °2  *  °*5 j  -  0. 


NAA 


p~  M  •  N\/NT(  p)  p^  Pm,n2  '  p=  0,1 ’2 . pmax  ’  (l3) 


where  NT(p)  is  the  number  of  power  terms  in  annular  band  p , 

S  n2  is  the  sum  of  the  squares  of  the  powers  in  annular 
band  p,  and  ^  is  the  normalization  factor.  (This  factor  is 
necessary  only  when  the  transform  of  a  function  F  has  been 
computed  by  the  FFT,  as  in  this  study.  This  is  because  the 
coefficients  from  an  FFT  appear  multiplied  by  M*N,  that  is, 
they  appear  unnormalized.  Whenever  equation  (6)  is  used  to 
compute  a  discrete  transform,  the  coefficients  are  already 
normalized  by  j^,  thus  eliminating  the  need  for  normalization 
by  (13).)  Therefore,  an  array  of  P  y+l  normalized  average 
amplitudes  NAAp  can  be  calculated  from  the  discrete  Fourier 
transform  of  a  cloud  sample.  When  these  NAA's  are  arranged 
in  order  of  increasing  p  ,  the  resulting  vector  is  called  a 
power  spectrum. 

For  a  3?  x  37  transform  array  the  annular  integration 
technique  reduces  1,369  Fourier  coefficients  to  a 


r  1 

P  +  1  =  / 

max  L  | 

spectrum  (a  re 
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In  order  to  better  understand  the  concept  of  annular 
integration  of  power,  consider  this  example.  Let  M  =  N  =  8. 
The  number  of  distinct,  concentric  annular  bands  generated 

by  (10)  is  r  t=r ~ - r~T2~  1  r  1 

V  f  +  [fj  +  0 . 5 J  +  l  0.5 J  +  1  =  7. 

Figure  6  shows  the  seven  annular  bands  of  an  8  x  8  discrete 
Fourier  transform. 

NAA  can  be  thought  of  as  a  measure  of  the  amplitude  of 

P 

the  harmonics  in  F  with  wavenumber  p  in  one  direction.  All 
coefficients  whose  powers  are  used  to  calculate  NAAp  contain 
at  least  one  subscript  m  or  n  which  is  equal  to  or  is  "nearly 
equal  to"  p.  Each  NAA  is  a  weighted  average  of  the  ampli¬ 
tudes  of  any  two-dimensional  waveforms  of  wavelength  either 

TUT 

"approximately"  —  pixels  in  the  x-direction  or  "approximately" 
N  p  ^ 

—  pixels  in  the  y-direction,  given  that  F  is  defined  for  all 
points  (i,j)  on  an  M  x  N  grid. 

For  a  37  x  37  array  of  satellite  data  values  which  repre¬ 
sent  a  25  x  25  nm  region  of  the  earth,  the  wavelength  of  a 

/\ 

harmonic  in  an  imagery  field  F  with  wavenumber  P  in  either  the 
x-direction  or  the  y-direction  is  inversely  proportional  to 
the  wavenumber: 

\  = 

pixels/cycle]  •  ^  nm/pixel^ 

-  p  nm/cycle.  (14) 

For  this  study,  therefore,  NAAp  is  the  average  amplitude  of 
all  two-dimensional  waves  (cloud  patterns)  whose  wavelengths 


=  pixels/cycle 

=  (Cpl 


25 

are  —g  nm 


in  either  the  x-direction  or  the  y-direction  of  the 

A. 

cloud  imagery  sample  F .  Figure  7  shows  wavelengths  as  a  func¬ 
tion  of  P. 

As  an  example,  it  is  easy  to  compute  the  power  spectra 
of  G  (i, j)  (7)  and  of  H  (i,j)  (8): 


6 

H 

NAAn  =  0 

NAA-  = 

0 

NAAj  =  20 

,  0 

NAAj  c 

0 

NAA2  =  0 

naa2  = 

0 

NAA_  =  0 

NAA-  = 

i 

3 

3 

4 

Note  that  only  one  NA^,  is  nonzero  in  each  case.  There  is 
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Figure  6.  The  seven  annular  bands  of  the  transform  of  an  8  x  8 
array  as  generated  by  (10).  Each  "block"  on  the  m-n  plane  repre 
sents  the  power  P  n  of  a  Fourier  coefficient  $  (see  (9)). 


Figure  7.  Wavelength  as  a  function  of  annular  hand  number 


enough  of  a  difference  between  these  two  spectra  to  be  able 
to  tell  which  spectrum  is  G's  and  which  spectrum  is  ft's. 
Such  differences  will  be  used  to  determine  a  sample's  cloud 
type  given  its  spectrum. 

It  is  interesting  to  note  that  annular  band  0  contains 

the  power  of  just  one  coefficient,  n.  From  the  FFT,  n  n 

/s  U ' U  i5D  .  /v  .  u  • u 

is  the  unnormalized  average  of  F  over  *71.  Since  F  is  real¬ 
valued,  Im(^Q  Q)  =  0.  Letting  p  =  0,  from  (13)  we  have: 


normalization  factor 


=  mean  of  F.  (15) 

One  important  property  of  the  annular  integration  tech¬ 
nique  is  that  it  is  not  sensitive  to  directionality.  Any 
finite  discrete  periodic  function  F(i,j)  on  rotated  any 
integer  multiple  of  ^  radians^  will  have  the  same  power  spec¬ 
trum  it  "originally"  had,  even  though  most  of  the  Fourier 
coefficients  change.  The  coefficients  in  a  Fourier  trans¬ 
form  are  sensitive  to  the  absolute  directionality  (orienta¬ 
tion)  of  the  harmonics  in  a  function,  whereas  a  power  spectrum 
is  not.  In  other  words,  a  power  spectrum  retains  information 
relative  to  the  shapes  and  sizes  of  clouds  in  an  imagery 

/N 

sample  F  but  loses  the  sense  of  their  directionality.  This 
is  an  advantageous  property  for  a  power  spectrum  to  have, 
since  the  north-south,  east-west  alignment  of  clouds  should 
not  affect  their  classification.  It  is  for  this  reason  that 
the  spectral  classifier  uses  the  imagery  sample's  power  spec¬ 
trum  to  identify  its  cloud  type. 

^One  need  only  be  concerned  with  90°  rotations,  since  any  other 
rotation  of  F  on  ould  result  in  F's  not  being  ^defined  for 
all  grid  points  in  &l.  Remember  F  may  change  but  ^remains 

fixed . 

25 


Figure  8  is  an  example  of  the  rotation  property.  Func¬ 
tion  (a)  is  rotated  90°  about  the  point  (1,1)  to  obtain  func 
tion  (b).  The  average  value  tji  q  q  of  each  function  is  2, 
but  in  (a)  the  only  other  nonzero  coefficient  is  t{)  ^  q  =  1, 
whereas  in  (b)  the  only  other  nonzero  coefficient  is 
4>0  4  =  1  •  It  is  easily  seen  that  ip  ^  q  (from  (a))  and 
1 1>0  ^  (from  (b))  are  used  to  calculate  each  function's  re¬ 
spective  NAA^,  since  (from  (10)): 


A2  +  o2  +  0.5  =  /o2  + 


4*  +  0.5J  =  4. 

The  patterns  in  functions  8(a)  and  8(b)  are  the  same  except 
for  their  orientation.  Their  power  spectra  are  equal,  even 
though  the  functions  and  their  transforms  are  different. 

An  infinite  number  of  test  patterns  could  have  the 
power  spectrum  given  by  Figure  8.  In  addition  to  the  cases 
given,  there  could  be  variations  in  both  the  i  and  j  direc¬ 
tions  since  any  of  the  22  Fourier  coefficients  in  p  =  4 
could  be  nonzero  so  long  as  the  sum  of  the  squares  of  their 
powers  is  equal  to  1 . 

Annular  integration  is  an  effective  data  reduction  tech 
nique  and  a  useful  tool  for  objectively  extracting  cloud 


size  information  from  an  imagery  sample  F.  As  with  any  data 
reduction  technique,  some  information  is  lost.  This  tech¬ 
nique  loses  the  sense  of  direction  of  the  cloud  elements. 

It  can  also  confuse  highly  directional  cloud  features,  e.g., 
cloud  bands,  with  more  complex  scenes  having  variation  in 
all  directions,  e.g.,  cumulus  clouds.  The  losses  are,  how¬ 
ever,  adequately  compensated  by  the  convenience  and  speed  of 
working  with  a  one-dimensional  power  spectrum.  Using  AFGL's 
CDC  6600  computer,  the  fast  Fourier  transform  algorithm  cal¬ 
culates  the  discrete  Fourier  transform  of  a  37  x  37  input 
array  F  (1,369  numbers)  in  7/10  sec,  and  the  calculation  and 
annular  sorting  and  summing  of  the  powers  of  these  coeffi¬ 
cients  are  done  in  3/100  sec.  Such  computational  speed  can 
allow  for  quick  cloud  sample  identification. 
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Figure  8.  Illustration  of  the  rotation  property  of  the  annular 
integration  technique.  Array  (a)  is  rotated  90°  about  (1,1)  to 
obtain  (b).  Note  the  pattern  in  (a)  occurs  four  times  horizont¬ 
ally  (m  *=  4)  and  is  unchanged  vertically  (n  ■  0 ) *  thus  ^  q  =  1 . 

Just  the  opposite  is  true  for  (b)  (m  ■  0  and  n  •  4),  hence  ^ 

■  1.  Although  oriented  differently  each  of  these  two  arrays0’ 
(functions)  exhibits  the  same  pattern.  The  spectral  classifier 
would  identify  each  sample  array  as  the  same  "type"  since  their 
spectra  are  similar. 
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3.2  The  Discriminant  Function 

That  part  of  the  classifier  which  determines  from  what 
type  of  cloud  a  given  power  spectrum  has  been  generated  is 
called  a  Bayes  decision  rule .  Simply  stated,  the  cloud' type 
chosen  by  the  decision  rule  to  be  contained  in  a  satellite 

/X 

imagery  sample  F  is  the  one  that  is  most  likely  to  generate 
that  sample's  observed  power  spectrum.  In  order  to  accomplish 
this,  the  rule  assumes  some  knowledge  of  the  spectra  of  every 
cloud  type  it  is  allowed  to  choose. 

The  decision  rule  is  allowed  to  classify  a  sample  as 
one  of  nine  cloud  types:  stratocumulus  (Sc),  stratus  (St), 
cumulus  (Cu),  clear,  altocumulus  (Ac),  altostratus/nimbo- 
stratus  (As/Ns)  ,  cirrus  (Ci),  cirrostratus  (Cs),  and  cumulo¬ 
nimbus  (Cb).  For  convenience,  these  will  be  referred  to  as 
cloud  types  i  =  1,2,3,...,9. 

Let  the  finite  set  Q  =  . of  the  nine 

states  of  nature  ui  be  the  set  from  which  the  classifier 

chooses  a  cloud  type,  or  in  other  words,  the  state  of  a  cloud 
8 

sample.  Since  the  state  of  nature  is  a  chance,  unpredic¬ 
table  event,  w  is  called  a  random  variable . 

The  spectral  classifier  uses  the  following  discriminant 
function  to  help  the  decision  rule  decide  what  kind  of  cloud 
(which  w ^)  a  sample  contains: 

*!<*>  *  -  !<*-■<!>'  E‘<x-  «*!>  -  |  In  2- -  | 

+  In  P(w.),  i  =  1,2,3 . 9.  U6) 

x  9 

where  x  is  an  r-component  column  vector  ,  is  the  r-compo- 

nent  mean  vector  of  class  i,  t  denotes  transpose,  is  the 

r  x  r  class-conditional  covariance  matrix  for  class  i, 

is  the  inverse  of  , 

i 

7 

These  two  cloud  types  are  merged  into  one  category  in  order 
to  be  consistent  with  the  3DNEPH. 

Q 

E.g.,w-  means  the  imagery  sample  contains  cumulus  clouds. 

9  • 

For  this  study  it  is  an  observed  26-component  power  spectrum. 


is  the  determinant  of  ^  * 

^-"*i  i 
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and  P( u^)  is  the  a  priori  probability  that  i  occurs. 

The  matrices  p^  and  and  the  a  priori  probabilities 

P(oj^),  must  all  be  calculated  from  the  sample  set  of  images 
to  be  classified  prior  to  any  sample  classification  process 
which  uses  (16).  The  mean  power  spectrum  p^  of  cloud  class 
i  is  simply  the  average  of  all  the  power  spectra  generated 
from  samples  containing  cloud  type  i.  The  matrix  is  com¬ 
posed  of  the  covariances  a  -v  of  NAA'!'  and  NAA?;,  which  belong 

Jk  J 

to  cloud  type  i.  Along  the  main  diagonal  of  /  „•  where 

^  2 

j  =  k  the  covariances  become  variances  of  the 

kth  component  NAA^  of  the  spectra  that  were  generated  from 
samples  containing  cloud  type  i.  P(  at)  is  simply  the  rela¬ 
tive  frequency  of  occurrence  of  cloud  type  is 

N. 

P(Wi)  =  ,  (17) 

where  is  the  number  of  samples  in  the  data  set  to  be 
processed  that  contain  cloud  type  i,  and  T  is  the  total  num¬ 
ber  of  imagery  samples  in  the  data  set.  Note  that  0  <  P(u>i) 

9 

^  1,  and  that  ^  P(w^)  =  1. 
i=l 

Equation  (16),  except  for  the  term  In  P(  a^) ,  is  merely 
the  natural  log  of  the  multivariate  normal  density s 

«  - i  (x-  p)ty-1(x-  p) 

nfxl  =  - - 7-TT  pL  2  ^  J.  Mfi) 


p(x)  =  - r7?\ — 17?  eL  J’  (1( 

(27r)r/2y 1/2 

where  x  ,  ^  ,  and  p  are  as  previously  defined.  Hence  d^ 
is  a  measure  of  the  probability  of  a  cloud  sample's  being 
in  state  g-ven  that  sample's  power  spectrum  x.  The  i 


for  which  d^  is  largest  is  the  most  reasonable  choice  of 
cloud  type  and  one  that  will  be  correct  most  of  the  time. 
Such  a  choice  would  make  p(x)  and,  therefore,  the  a  posteri¬ 
ori  probability  P(«^  |  x),  the  probability  of  cloud  type  i 
given  x,  largest. 

Let  x  be  a  26-component  observed  power  spectrum  gener¬ 
ated  from  a  37  x  37-pixel  cloud  sample  (thus  in  (16), 
r  =  26).  The  term  ^  In  2  it  =  13  In  2»  is  merely  an  additive 


t 


constant  which  is  independent  of  i,  and  can  therefore  be 
ignored.  Statistical  independence  between  the  NAA’s  of  the 
power  spectra  of  the  same  cloud  type  is  assumed.  Thus: 


/  2 

_  |°k  *  J 

"  (0  .3 


for  j,k  =  0,1,2... . ,25,  (19) 


so  that  reduces  to  a  26  x  26  diagonal  matrix  of  variances 
of  NAA^,  k  =  0,1,2, ...,25.  Rewriting  (16)  with  the  above 
conditions  yields:  -  ^ 

di^  =  “  2  [Co  0 ^  ^X1  “  ^1)-“(x25  “  ^25^  ^  °**# 


( 45y 


xo  -  Vo 
X1  - 


:  .  ■  I ln  Si  +  ln  p(“i) 


X25  “  M25 


or,  simplifying, 

25  i  2 

d  i  (*)  =  -  4  X!  ■  !n~~  -  I  In  Si  +  In  P(»i), 

1  2  n=0  (o^)2  2 

i  =  1,2,3 . 9,  (20) 

where  x  is  an  observed  power  spectrum,  a  26-component  column 
vector  comprised  of  xQ,  Xj ,  x^>  •••»  *25  (xnr  NAAn)s  ^  *s 
the  mean  of  all  NAA  ,  class  i  (from  |».)s  (o^)^  is  ^e 
variance  of  NAA^,  class  i  (from  £^);  P(w^)  is  the  a  priori 
probability  (17).  A  d^  is  calculated  for  each  of  the  nine 
cloud  types  Sc,  St,  Cu ,  ...,  Cb.  The  decision  rule  decides 
class  k  if  d^  >  d^  for  all  i  /  k. 

In  the  form  (20)  it  is  easy  to  see  that  a  cloud  imagery 
sample's  observed  power  spectrum  x  is  compared  component-by¬ 
component  to  the  mean  power  spectra  of  each  cloud  type  by 

i  2 

measuring  the  squared  distance  (xn  -  nn)  between  xn  (NAAn) 


and  for  all  n,  all  i.10  The  smaller  (xn  -  n^)2  is,  the 
larger  d^  is,  so  that  the  i  for  which  d^  is  largest  indicates 
the  cloud  type  whose  mean  or  "characteristic"  power  spectrum 
most  closely  resembles  x . 

This  is  how  the  decision  rule  uses  the  discriminant 
function  (20)  to  make  its  classification  decision.  The  auto¬ 
mated  spectral  classifier  is  for  this  reason  called  a  minimum 
distance  classifier .  If  each  p^  is  thought  of  as  an  ideal 
typical  spectrum  for  cloud  class  i,  then  (20)  may  also  be 
described  as  a  template-matching  procedure. 

In  cases  where  x  is  "close"  to  more  than  one  character¬ 
istic  power  spectrum  p^  (the  left-most  term  on  the  right  side 
of  (20)  is  or  is  nearly  the  same  for  more  than  one  cloud 
type  i),  the  classifier's  final  decision  depends  on  P(uk). 

For  example,  when  x  is  equally  (or  nearly  equally)  close  to 
two  mean  spectra  p^  and  p^,  the  choice  made  by  the  classifier 
favors  the  a  priori  more  likely  cloud  class.  So,  if 
P(«k)  >  P(«i),  then  dk  >  d^  and  cloud  type  k  would  be  chosen. 

A  priori  probabilities  may  be  used  to  more  suitably 
adapt  the  classifier  to  different  regions  of  the  world.  For 
instance,  in  equatorial  regions  P(u^)  =  P(Cb)  is  much  greater 
than,  say,  P((it>)  =  P(St).  However,  in  maritime  midlatitude 
areas,  just  the  opposite  is  likely  the  case.  Differences  in 
relative  frequency  of  occurrence  of  cloud  types  over  various 
geographical  areas  can  be  reflected  in  a  proper  choice  of 
P(w  i) 's  for  each  of  those  areas.  This  is  an  advantageous 
versatility  of  the  spectral  classifier  technique. 

The  discriminant  function  can  be  thought  of  as  a  mathe¬ 
matical  formalization  of  common-sense  rules.  It  is  a  logi¬ 

cal,  effective,  and  straightforward  way  of  determining  the 
cloud  type  of  an  imagery  sample  given  that  sample's  power 
spectrum. 

In  addition  to  using  the  discriminant  function  (20)  to 
identify  each  cloud  sample,  two  other  variations  of  (20)  were 

*°The  numerator  of  the  left-most  term  on  the  right  side  of 
(20)  can  be  written  as  (  x  -  (x-  p^)  *  ||  x  -  p^||. 

This  is  called  the  Euclidean  distance  from  x  to  p^. 

V 


utilized  to  determine  how  different  parts  of  (20)  affect 
classification  decisions. 

The  mean  classifier,  which  bases  its  cloud  type  deci¬ 
sions  solely  on  the  mean  visible  and  infrared  values  (NAAQ) 
of  the  cloud  samples,  simply  computes  the  linear  distance  be¬ 
tween  two  points: 

m.(VIS  .  ,  IR  ,  )  =  /(vTs  .  -  vTs.)2  +  (IR  .  -  IR.)2, 

1  obs  obs  vv  obs  x  obs  1 

i  =  1,2,3 . 9,  (21) 

where  vlsol;)S  and  IR0bs  are  olDserved  average  visible  and 

infrared  values  of  the  sample  being  classified,  and  VIS^  and 

IR^  are  the  mean  visible  and  infrared  values  of  all  samples 

of  cloud  type  i.  Before  (21)  can  be  used,  VIS.  and  IR.  must 

11  .  1  1 

be  known  for  each  cloud  class  i.  An  m  is  calculated  for 
each  i,  and  cloud  type  k  is  chosen  if  m^  <  irn  for  all  i  /  k. 
The  results  obtained  by  this  classification  procedure  will 
indicate  how  much  the  means  of  a  sample  contribute  to  cloud 
type  classification. 

Another  variation  of  (20),  which  is  essentially  a  curve¬ 
fitting  operation  in  a  least  squares  sense,  bases  its  deci¬ 
sions  solely  on  the  means  and  variances  of  the  power  spectra 
of  each  cloud  typei 

j  n  (Jtn  ■  M*)2 

«.(*)  =  -  4  y*  — ■  --0  ■ —  ,  i  =  1,2,3 . 9,  (22) 

1  C  r,  /  1 

n=0  (0n; 

where  x ,  x^,  and  (o^)2  are  as  defined  for  (20).  Equation 

(22)  is  merely  equation  (20)  without  the  a  priori  probability 
and  covariance  terms.  As  in  (20),  the  observed  power  spectrum 
x  is  compared  component-by-component  to  each  mean  spectrum 
by  calculating  a  6  for  each  of  the  nine  cloud  classes. 

Cloud  type  k  is  chosen  if  6^  >  6^  for  all  i  /  k.  When  com¬ 
pared  with  the  results  of  (20),  the  performance  rate  of  (22) 
will  show  how  much  the  classifier  relies  on  a  priori  pro¬ 
babilities  . 


The  mean  classifier  follows  the  same  logic  as  does  the 
spectral  classifier. 


3-3  Data  Availability  and  Classifier  Versatility 


During  nighttime  orbits  only  infrared  cloud  data  are  ob 
tained  by  satellites.  There  are  also  times,  albeit  less  fre 
quent,  when  only  visible  data  are  available  for  use.  Never¬ 
theless,  more  often  than  not  both  visible  and  infrared  data 
are  available  for  use  during  the  daytime  scans  of  a  polar 
orbiting  satellite. 

In  order  to  adapt  the  classifiers  (20),  (21),  and  (22) 
to  each  of  the  above  conditions,  simple  and  straightforward 
modifications  were  made  to  each. 

The  visible-only  and  infrared-only  adaptations  of  the 
mean  classifier  (21)  are  defined  as,  respectively, 

mYIS  »  |vTSobs  -  VTS.|  ,  i  =  1,2,3 . 9,  (21a) 

and  mi:[-R  =  \™ohs  -  IR^--*  i  =  1,2,3 . 9,  (21b) 

where  VISobs>  VIS^,  IRobs ,  and  IR^  are  as  defined  for  (21). 
Formulas  (21a)  and  (21 b)  make  a  cloud  type  decision  the  same 
way  as  (21)  does. 

The  discriminant  functions  d^(x)  (20)  and  6^(x)  (22) 
were  adapted  to  separate  and  simultaneous  visible/infrared 
processing  in  the  same  way.  The  spectral  classifiers  (20) 
and  (22)  use  a  visible  power  spectrum  x,  generated  from  a 
visible  imagery  sample,  as  input  for  the  visible-only  condi¬ 
tion,  and  use  an  infrared  spectrum  y  for  the  IR-only  condi¬ 
tion.  Using  (22)  as  an  example,  when  both  visible  and 
infrared  imagery  are  available  simultaneously  with  which  to 
classify  a  cloud  sample,  the  individual  terms  5 ^(x) 

(visible)  and  6  ^(  y)  (infrared)  are  calculated  separately 
for  each  i  and  then  added  together.  The  cloud  type  i  for 
which  the  sum  &  ^  (x)  ♦  i.(y)  is  largest  is  chosen. 
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4 .  RESULTS 


The  available  DMSP  5D  visible  and  infrared  imagery 
was  displayed  on  a  McIDAS  terminal  TV  screen.  Cloud  imagery 
samples  were  then  interactively  found,  examined,  and  sub¬ 
jectively  classified  by  skilled  analysts.  Samples  contain¬ 
ing  mixed  clouds  were  avoided,  as  were  those  with  coastal 
backgrounds  and  snow-covered  areas.  Once  collected  on  tape, 
a  group  of  cloud  samples  would  then  be  classified  objectively 
using  the  spectral  classifier  program. 

The  classifier  was  allowed  to  run  on  two  separate  sets 
of  data.  The  first  set,  which  contained  73  cases,  was  sub¬ 
jectively  classified  by  a  different  analyst  than  was  the 
second  set  of  70  samples. 

The  overall  performance  of  the  classifier  on  both  sample 
sets  was  highly  similar  (see  Table  1).  This  suggests  that 
the  classification  technique  is  sound  and  consistent.  For 
these  reasons,  the  results  of  both  trials  have  been  combined 
in  the  following  discussion. 

TABLE  1.  CLASSIFICATION  ACCURACY  FOR  TWO  DATA  SETS 


Visible 

IR 

Combined 
Vis  and 

Mean  Classifier, 

Data  Set  1 
(73  samples) 

34% 

37% 

49% 

Mean  Classifier, 

Data  Set  2 
(70  samples) 

37% 

43% 

54% 

Spectral  Classifier, 
Subset  1 

66%> 

68% 

82% 

Spectral  Classifier 
Subset  2 

64% 

63% 

79% 

The  composition  of  the  combined  set  of  cloud  samples  is 
shown  in  Figure  9-  Cirrus  is  the  most  common  cloud  type  in 
the  sample  set  with  36  of  the  143  cases.  Cumulus  is  the  next 
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TOTAL  SAMPLES- 143 


Figure  9.  Composition  of  the  cloud  imagery  sample  set 


most  frequent  with  24  samples.  Mid-level  clouds  (Ac,  As/Ns) 
did  not  appear  too  often  within  the  available  imagery  and  thus 
have  a  relatively  low  number  of  cases:  9  for  Ac  and  5  for 
As/Ns.  The  backgrounds  of  these  cloud  samples  range  from  the 
Mediterranean  Sea  to  the  Balkan  peninsula,  and  the  deserts  of 
the  Middle  East  and  northeast  Africa. 

With  the  information  shown  in  Figure  9.  the  a  priori 
probabilities  P(u>^)  (17)  can  be  computed  for  each  of  the  nine 
cloud  types.  If  an  observer  had  to  classify  this  sample  set 
on  the  basis  of  this  information  only,  then  optimum  results 
would  be  achieved  by  choosing  the  most  frequent  cloud  type  in 
the  sample  set. 

In  this  case,  "blindly"  choosing  cirrus  as  the  cloud 
type  of  all  143  imagery  samples  would  result  in  a  25$  correct 
classification  rate  (P(o)^)  =  P(Ci)  =  .251),  since  the  analyst 
is  sure  that  36  of  the  143  samples  are  cirrus.  This  is  con¬ 
siderably  higher  than  a  totally  random  classification  approach 
which  yields,  on  the  average,  p  =  1/9  =  .111,  or  an  11$  cor¬ 
rect  classification  rate. 

The  overall  results  of  the  three  classifiers  (20),  (21), 
and  (22)  should  be  compared  to  the  above  "blind"  performance 
rates  as  a  measure  of  the  classifiers’  effectiveness  in  deter¬ 
mining  cloud  type. 

4.1  Mean  Classifier 

The  performance  rate  of  the  mean  classifier  was  better 
than  that  of  the  "blind"  classification  approach  (see  Table 
2a).  Attempting  to  identify  from  one  of  the  nine  cloud  types, 
the  mean  classifier  classified  36$  of  the  cloud  samples  cor¬ 
rectly  using  visible  data  only  (2la),  40$  of  the  cloud  samples 
using  infrared  only  (21b),  and  52$  of  the  cloud  samples  using 
visible  and  infrared  data  (21).  There  was  a  rather  large  num¬ 
ber  of  misclassif ications  when  this  method  was  utilized, 
which,  although  an  improvement  over  the  25$  of  the  "blind 
observer,"  implies  that  the  mean  of  a  cloud  sample  itself  is 
not  sufficiently  discriminatory  to  clearly  distinguish  nine 
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TABLE  2a.  MEAN  CLASSIFIER  ACCURACY  FOR  NINE  CATEGORIES 


Cloud  Type 

Number  of 
Samples 

Visible 

IR 

Combined 
Vis  and 

Sc 

7 

29$ 

29$ 

29$ 

St 

10 

50$ 

20$ 

60$ 

Cu 

24 

42$ 

58$ 

58$ 

Clear 

19 

k2% 

74$ 

58$ 

Ac 

9 

11% 

0 

44$ 

As/Ns 

5 

60% 

40$ 

60$ 

Ci 

36 

25$ 

17$ 

42$ 

Cs 

19 

47$ 

42$ 

63$ 

Cb 

14 

29$ 

64$ 

50$ 

Total 

143 

36$ 

40$ 

52$ 

TABLE  2b.  MEAN 

CLASSIFIER 

ACCURACY  FOR 

FIVE 

CATEGORIES 

Cloud  Group 

Number  of 
Samples 

Visible 

IR 

Combined 
Vis  and 

Low 

41 

54$ 

66$ 

71$ 

(Sc/St/Cu) 

Mid 

14 

36$ 

36$ 

57$ 

(Ac/As/Ns) 

High 

55 

44$ 

25$ 

53$ 

(Ci/Cs) 

Clear 

19 

42$ 

74$ 

58$ 

Cb 

14 

29$ 

64$ 

50$ 

Total 

143 

44$ 

48$ 

59$ 

cloud  types  in  a  consistent  manner.  Figure  10  has  a  plot  of 
visible  and  infrared  averages  for  each  cloud  type.  Some  of 
the  averages  are  so  close  (Ac,  As/Ns,  St)  that  misclassif ica- 
tions  are  expected. 

A  significant  factor  reflected  in  these  results  is  the 
diverse  range  of  backgrounds  found  among  the  samples.  This 
greatly  affects  the  means  of  samples  of  the  same  cloud  type, 
especially  those  in  which  a  lot  of  background  shows  through 
(Cu,  clear,  Ci).  For  example,  one  Cu  case  has  a  mean  visible 
count  of  19.  while  another  has  one  of  109,  almost  six  times 
as  high. 

When  the  cloud  samples  are  classified  into  fewer  cate¬ 
gories,  overall  results  improve.  The  nine  cloud  types  were 

grouped  into  five:  low  (Sc,  St,  Cu),  middle  (Ac,  As/Ns), 

1 2 

high  (Ci,  Cs),  clear,  and  Cb.  Thus,  for  example,  a  Sc 
sample  classified  as  Cu  was  considered  a  correct  identifica¬ 
tion,  or  "hit." 

Using  visible  means  only,  low  clouds  were  correctly 
hit  54%  of  the  time,  middle  clouds  only  36%  of  the  time,  and 
high  clouds  44%  of  the  time.  For  all  five  categories,  44% 
were  classified  correctly. 

Using  infrared  only,  48%  hits  were  achieved:  66%  low 
clouds  were  hit,  36%  middle  clouds  were  hit,  but  (surpris¬ 
ingly)  only  25%  high  clouds  were  correctly  classified. 

When  visible  and  infrared  means  were  used  simultane¬ 
ously,  71%,  57%.  and  53%  of  the  low,  middle,  and  high  cloud 
samples,  respectively,  were  correctly  classified  for  a  com¬ 
bined  five-category  59%  hits. 

If  the  above  five-category  results,  listed  in  Table  2b, 
are  compared  with  the  "blind"  classification  approach  per¬ 
formance  rate  of  38%  (attainable  by  calling  each  sample 

1 2 

Cs  and  Cb  were  confused  aV  times  by  the  classifiers  since 
they  often  look  very  similar  in  the  visible  and  in  the  infra¬ 
red.  Because  of  the  severe  weather  frequently  associated  with 
Cb,  it  was  felt  that  these  two  cloud  types  should  be  in  sepa¬ 
rate  classes  even  though  results  would  improve  if  they  were 
grouped  together. 
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Figure  10.  Average  visible  and  infrared  means  for  each 
cloud  type,  sample  subset  2  (70  samples).  High  values  on 
the  infrared  scale  represent  cold  temperatures. 


"high"),  then  it  is  seen  that  the  mean  classification  approach 
has  an  advantage  over  chance. 

4.2  Average  Power  Spectra  Characteristics 

This  section  discusses  the  results  obtained  when  the 
cloud  type  decision  rules  (20)  and  (22)  were  utilized  on  the 
samples.  Before  these  two  methods  could  be  applied,  the  means 
and  variances  of  the  power  spectra  of  each  cloud  type  were 
computed . 

The  average  visible  and  infrared  power  spectra  of  each 
of  the  nine  cloud  types  contained  in  sample  set  2  ( 70  cases) 
are  plotted  in  Figures  11  and  12.  As  can  be  seen  in  these 
graphs,  the  most  noticeable  differences  among  the  average 
spectra  lie  in  the  lower  wavenumbers,  and  especially  in  the 
means  (NAAQ  (15))  of  the  cloud  types.  Cs  has  the  highest 
visible  and  infrared  means,  while  Cb  has  the  second  highest. 
Appropriately,  clear  has  both  the  lowest  visible  and  lowest 
infrared  means.  (High  values  on  the  infrared  scale  corres¬ 
pond  to  low  temperatures  while  low  values  correspond  to  high 
temperatures. ) 

As  previously  stated,  coefficients  of  a  Fourier  expan¬ 
sion  are  the  amplitudes  of  waves  which  lie  about  the  mean  of 
the  function  they  describe  (see  Figure  3)»  They  are  therefore 
some  measure  of  how  much  the  values  of  that  function  deviate 

from  -its  mean.  Hence  coefficients  i  „  that  lie  in  the 

m ,  n 

smaller-radius  ( lower  p )  annular  bands  are  the  amplitudes 
of  larger-scale  variations  (larger  X  -  see  equation  (14)) 
in  cloud  brightness  and  spatial  configuration,  and  the  coeffi¬ 
cients  that  lie  in  the  higher  p  annular  bands  are  the  ampli¬ 
tudes  of  smaller-scale  variations  (smaller!  ).  It  is  in 
this  sense  that  the  average  power  spectra,  generated  from 
Fourier  coefficients,  are  indicative  of  the  amount  of  texture 
present  in  the  imagery  field.  As  the  higher  wavenumber  NAA's 
increase,  the  imagery  sample  appears  less  uniform. 

Beyond  wavenumber  6,  the  range  of  values  (the  differ¬ 
ence  _  NAAm^n)  for  the  average  spectra  of  the  strati- 
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Figure  12.  Semi-log  plot  of  infrared  average  power  spectra 
of  each  cloud  type,  sample  subset  2. 


form  cloud  types  (As/Ns,  St,  Cs)  is  particularly  small, 
especially  in  the  infrared,  indicating  that  their  imagery 
is  quite  smooth  in  appearance.  On  the  other  hand,  the  corres¬ 
ponding  range  of  the  cumuliform  cloud  type  spectra  (Cu,  Sc, 

Ac)  is  significantly  larger,  especially  in  the  visible.  J 
This  implies  that  cumuliform  imagery  is  much  less  uniform 
in  appearance  than  is  stratiform  imagery.  Figures  13  and 
14  are  plots  of  average  spectra  for  wavenumbers  greater  than 
6. 

4.3  The  Discriminant  Functions 

The  Bayesian  decision  rules  (22)  and  (20),  as  expected, 
performed  much  better  than  the  mean  classifier.  Improve¬ 
ments  of  20  to  30$  over  the  mean  classifier  performance  rates 
were  obtained  when  the  discriminant  functions  were  used. 

Such  a  significant  increase  in  classifier  accuracy  is  mostly 
due  to  the  fact  that  there  is  a  lot  more  cloud  feature  infor¬ 
mation  present  in  the  power  spectrum  of  a  cloud  sample  than 
there  is  in  just  its  mean.  Results  are  shown  in  Tables  3 
and  4. 

The  largest  increase  observed  was  one  of  32$.  Using 
visible  data  only,  the  accuracy  of  the  mean  classifier  was 
36$,  whereas  that  of  the  minimum  distance  classifier  (20) 
was  68$. 

In  general,  Cs  was  the  cloud  type  most  often  correctly 
classified,  especially  by  (20).  When  both  visible  and  infra¬ 
red  data  were  used  to  determine  cloud  type,  1 7  of  19  Cs 
cases  were  properly  identified  (89$) •  Cb  cloud  samples  were 
also  well  classified  by  both  (20)  and  (22). 

Mid-level  cloud  types  were  the  most  inconsistently 
classified.  For  example,  only  two  Ac  samples  were  called  Ac 
by  (22)  when  visible  was  used  alone.  However,  when  visible 

^Stratiform  range  p  >  6?  Vis?  «  .40,  IR s  «  .04. 

Cumuliform  range  p  >6j  Vis?  «  1.0,  IR;  «  .40 
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Figure  13-  Linear  plot  of  visible  average  power  spectra  of 
each  cloud  type  beyond  annular  band  6,  sample  subset  2. 
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Figure  14.  Linear  plot  of  infrared  average  power  spectra  of 
each  cloud  type  teyond  annular  band  6,  sample  subset  2. 


TABLE  3.  NINE-CATEGORY  SPECTRAL  CLASSIFIER  ACCURACY  WITHOUT 
A  PRIORI  AND  COVARIANCE  TERMS  (22) 


Cloud  Type 

Number  of 
Samples 

Visible 

IR 

Combined 
Vis  and  IR 

Sc 

7 

57# 

29# 

86# 

St 

10 

40# 

20# 

70# 

Cu 

24 

67% 

54# 

79# 

Clear 

19 

7k% 

74# 

74# 

Ac 

9 

22% 

33# 

67# 

As/Ns 

5 

20% 

0 

100# 

Ci 

36 

42# 

72# 

83# 

Cs 

19 

74# 

84# 

89# 

Cb 

14 

79# 

79# 

86# 

Total 

143 

57% 

61# 

81# 

TABLE  4.  NINE- 

-CATEGORY  SPECTRAL  CLASSIFIER  ACCURACY  WITH 

A  PRIORI  AND  COVARIANCE  TERMS  (20) 

Cloud  Type 

Number  of 
Samples 

Visible 

IR 

Combined 
Vis  and  IR 

Sc 

7 

43# 

43# 

71# 

St 

10 

70# 

80# 

100# 

Cu 

24 

58# 

50# 

71# 

Clear 

19 

95% 

89# 

89# 

Ac 

9 

78# 

44# 

78# 

As/Ns 

5 

100# 

100# 

100# 

Ci 

36 

44# 

72# 

64# 

Cs 

19 

84# 

74# 

89# 

Cb 

14 

79# 

71# 

86# 

Total 

143 

68# 

69# 

79# 
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alone  was  used  by  (20),  7  of  9  Ac  cases  were  correctly 
classified.  A  similar  statement  can  be  made  of  As/Ns  sam¬ 
ples  using  their  infrared  spectra  only  (0  out  of  5  vs.  5  out 
of  5).  Of  course,  performance  of  the  classifier  on  Ac  and 
As/Nd  cloud  types  should  not  be  judged  as  either  poor  or  good 
on  the  basis  of  this  information  alone,  since  only  a  small 
number  of  these  mid-level  samples  were  selected  and  processed. 

On  the  other  hand,  clear  samples  were  most  consistently 
classified  by  the  discriminant  functions,  as  can  be  seen  in 
Table  3  (74#)  and  Table  A  (89#). 

The  mistakes  the  classifier  made  when  it  erred  were  of¬ 
ten  reasonable.  Oftentimes  Sc  was  mistaken  for  Cu,  or  Cs 
for  Cb.  A  measure  of  how  often  such  mistakes  occur  lies  in 
comparing  the  discriminant  function  nine-category  results 
with  the  discriminant  function  five-category  (low,  middle, 
high,  clear,  Cb)  results. 

The  overall  five-category  combined  performance  rates 
increased  up  to  10#  over  those  of  the  nine-category  set. 

This  implies  that  a  substantial  number  of  samples  misclassi- 
fied  as  one  of  nine  types  were  not  misclassified  as  one  of 
five  types.  Results  are  shown  in  Tables  5  and  6. 

The  improvement  of  the  accuracy  of  the  discriminant 
functions  when  allowed  to  classify  from  five  types  instead 
of  nine  should  not  be  misinterpreted.  Five-category  accu¬ 
racies  cannot  get  worse.  By  retabulating  the  results  for 
nine  cloud  categories  into  five,  the  overall  results  of  the 
classifier  can  only  remain  the  same  or  get  better.  There¬ 
fore,  the  important  point  is  not  whether  the  classifier  has 
improved,  but  by  how  much  it  did  improve.  Such  a  comparison 
indicates  the  amount  of  "reasonable"  errors  the  classifier 
makes . 

The  performance  of  the  visible-only  nine-category 
classifier  (22),  poorest  of  all  at  57$  (see  Table  3)*  im¬ 
proved  by  nearly  10#  to  66#.  This  was  better  than  any  of 
the  other  observed  increases. 

The  classifier  (20)  called  a  sample  from  the  low  cloud 
category  either  Sc,  St,  or  Cu  80#  of  the  time,  using  infra- 


TABLE  5.  FIVE-CATEGORY  SPECTRAL  CLASSIFIER  ACCURACY  WITHOUT 


A 

PRIORI  AND 

COVARIANCE 

TERMS 

(22) 

Cloud  Group 

Number  of 
Samples 

Visible 

IR 

Combined 
Vis  and  IR 

Low 

(Sc/St/Cu) 

41 

71% 

54$ 

83$ 

Mid 

(Ac/As/Ns ) 

14 

21% 

21$ 

79$ 

High 

(Ci/Cs) 

55 

69% 

84$ 

89$ 

Clear 

19 

74$ 

74$ 

74$ 

Cb 

14 

79% 

79% 

86$ 

Total 

143 

66% 

67% 

84$ 

TABLE  6.  FIVE-CATEGORY  SPECTRAL  CLASSIFIER  ACCURACY  WITH 

A 

PRIORI  AND 

COVARIANCE 

TERMS 

(20) 

Cloud  Group 

Number  of 
Samples 

Visible 

IR 

Combined 
Vis  and  IR 

Low 

(Sc/St/Cu) 

41 

66% 

80$ 

85$ 

Mid 

(Ac/As/Ns ) 

14 

86% 

64$ 

86$ 

High 

(Ci/Cs) 

55 

75% 

75$ 

85$ 

Clear 

19 

95% 

89$ 

89$ 

Cb 

14 

79% 

71$ 

86$ 

Total 

143 

76% 

77$ 

86$ 
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red  spectra  .  It  also  correctly  classified  three  of  ?  Sc  sam¬ 
ples,  8  of  10  St  samples,  and  12  of  24  Cu  samples  from  nine 
cloud  types  using  infrared  spectra.  Of  the  18  samples 
(4  Sc,  2  St,  and  12  Cu)  that  were  not  correctly  classified 
from  nine  types,  10  of  them  were  placed  within  the  low  cloud 
group.  Thus  more  than  half  of  the  previous  incorrect  classi¬ 
fications  were  recovered  because  of  five-category  partitioning. 

Although  high  clouds  had  the  highest  hit  percentage,  in 
general  both  high  and  low  clouds  had  similar  performance 
rates.  As  was  the  case  with  the  nine  cloud  types,  mid-level 
clouds  were  inconsistently  classified. 

4.4  Significance  of  A  Priori  Knowledge  and  Its  Effect  on  the 
Classifier 


Equation  (20)  differs  from  (22)  by  the  "a  priori"  terms 

-I  ln|£i|  and  In  P(w^).  These  extra  terms  are  intended  to 

discriminate  between  two  or  more  cloud  types  whose  mean 
spectra  are  all  "close"  to  an  observed  sample  spectrum  x. 
The  difference  in  the  results  of  using  (20)  and  (22)  was 
most  significant  when  using  visible  data  only.  Increases  of 
11  and  10$  were  observed  for  the  nine  category  (57  to  68%)  and 
five  category  (66  to  76$)  sets.  All  other  percentages  in¬ 
creased  on  the  order  of  2-8$. 

The  degree  of  help  that  the  a  priori  terms  give  the 
classifier  should  not  be  taken  for  granted  since  it  may  be 
unique  to  this  data  sample  set.  More  than  half  of  the  two 
largest  increases  of  11  and  10$  were  due  to  a  64$  increase  in 
the  number  of  mid-level  clouds  correctly  classified  when  (20) 
was  used.  As  described  earlier,  mid-level  clouds  were  the 
most  inconsistently  classified  group,  so  that  such  a  marked 
increase  should  be  carefully  considered  before  determining 
the  significance  of  these  results. 

In  those  cases  when  the  classification  rate  of  the  mid¬ 
level  clouds  did  not  change  so  abruptly  (IR-only,  combined 
visible  and  IR),  the  more  modest  improvements  of  2-8$  do  not 


overwhelmingly  suggest  the  use  of  (20)  over  (22).  It  is  true 
that  the  results  obtained  by  (20)  were  almost  always  better, 
but  the  processing  of  a  larger  data  set  would  provide  better 
information  with  which  to  determine  the  most  favorable  dis¬ 
criminant  function;  The  a  priori  probability  and  covariance 
terms  in  (20)  for  such  a  data  set  would  tend  to  be  a  better 
climatology  of  cloud  types. 

4.5  The  Classifier  vs.  Chance 

The  spectral  classifier  is  clearly  advantageous  over  a 
chance  classification  approach.  Choosing  the  highest  a  pri¬ 
ori  probability  from  nine  cloud  types  yields  a  25$  correct 
classification  rate  while  the  automatic  classifier  yields  an 
81$  correct  classification  rate.  Similarly,  choosing  the 
highest  a  priori  probability  from  five  cloud  types  yields  a 
38$  correct  classification  rate  while  the  automated  classifier 
yields  an  86$  correct  classification  rate. 

In  a  previous  similar  study,  Blackman  and  Pickett  (1 979) 

obtained  a  46$  classification  rate  for  a  six-category  set 

consisting  of  the  Sc,  Ac,  Ci,  Cs,  and  Cb  cloud  types.  Only 

DMSP  visible  data  of  1/3  ran  resolution  (along  the  satellite 
14. 

subtrack  )  were  available  for  this  study,  and  only  the  dis¬ 
criminant.  function  (20)  was  used.  Considering  only  these  six 
cloud  types,  this  study's  corresponding  correct  classification 
rate  (using  2/3  ran  visible  data  only,  along  with  (20))  was 
62$. 


4.6  Another  Data  Reduction  Technique 

The  question  of  how  little  information  is  actually 
necessary  to  enter  into  a  decision-making  process  in  order 
to  achieve  desirable  results  is  highly  important  since  prac¬ 
tical  considerations  require  that  the  great  volume  of  satel- 

14 

Refer  to  Bunting  and  Fournier  (I98O)  for  a  more  complete 
description  of  DMSP  satellite  sensor  capabilities. 
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lite  data  be  substantially  reduced.  This  technique  has 
helped  that  problem  by  reducing  the  number  of  data  values 
from  1,369  (an  input  imagery  sample)  to  26  (its  power  spec¬ 
trum).  However,  there  are  variations  of  this  method  that  can 
further  reduce  the  amount  of  information  it  needs  to  effec¬ 
tively  classify  cloud  samples. 

It  would  be  unwise  to  discard  any  parts  of  the  original 
imagery  array  itself,  since  the  transform  of  a  space-depen¬ 
dent  cloud  image  with  "holes”  in  it^  could  not  be  computed. 

But  once  the  imagery  is  transformed  into  its  Fourier  coeffi¬ 
cients,  the  following  two  possibilities  are  left  open:  (1) 
utilize  less  than  26  components  of  a  power  spectrum  in  the 
final  decision  process  (discriminant  function)?  and/or  (2) 
decrease  the  number  of  coefficients  used  in  computing  a 
cloud  sample's  power  spectrum. 

The  first  option  requires  a  more  detailed  study  of  the 
power  spectra  generated  from  the  transforms  of  the  imagery 
samples  of  each  cloud  type.  Blackman  and  Pickett  (1979) 
have  suggested  that  most  of  the  class-discriminating  power 
lies  in  only  a  few  of  the  26  spectral  components,  so  that  not 
all  of  them  are  necessary  to  obtain  the  results  already  pre¬ 
sented.  This  possibility  was  not  studied  but  it  deserves 
added  attention. 

The  second  option  was  studied  by  arbitrarily  decreas¬ 
ing  the  number  of  coefficients  in  computing  a  power  spectrum. 

Only  those  coefficients  i  with  non-negative  subscripts 

in  y  11 

m,n  >  0  were  allowed  into  any  power  spectra  calculatior s . 

All  others  were  ignored.  (A  more  complete  description  is 
given  in  Appendix  C.)  This  caused  minimal  discriminant 
function  modifications,  since  a  full  26-component  spectrum 
could  still  be  generated  from  these  fewer  coefficients.  By 
restricting  the  coefficients  used  to  those  which  lie  in  the 
first  quadrant  of  the  principal  part  of  the  wavenumber  plane 
(m,n  >.  0?  see  Figure  6),  at  least  one  Fourier  coefficient 
lies  in  any  given  annular  band  generated  by  (10).  This  im¬ 
plies  that  at  least  some  information  about  every  harmonic 

( i ,  j )  points  where  the  input  function  F(i,j)  would  be  undefined. 


(wavenumber)  present  in  the  cloud  imagery  field  will  be  in¬ 
cluded  in  its  spectrum. 

This  reduction  method  cut  the  number  of  Fourier  coeffi¬ 
cients  processed  for  the  power  spectrum  of  a  37  x  37  array 
from  1,369  to  36I,  a  drop  of  almost  4:1.  The  advantages  here 
are  computational  speed  and  storage  reduction,  but  the  dis¬ 
advantage  may  be  loss  of  information.  The  savings  in  compu¬ 
tational  time  are  small,  since  annular  integration  is  much 
quicker  than  computing  the  FFT.  Thus  a  slight  to  no  de¬ 
crease  in  the  performance  rate  of  this  technique  with  respect 
to  the  "all-coefficient"  method  would  justify  its  use. 

Results  for  the  first  quadrant  coefficient  method  are 
shown  for  both  the  nine  and  the  five  category  trials  in 
Tables  7.  8,  9.  and  10.  These  Tables  are  cross-referenced 
to  their  corresponding  "all-coefficient"  method  results.  The 
number  in  parentheses  which  immediately  follows  each  trial’s 
overall  nine-  or  five-category  performance  rate  (e.g.,  60% 
{+3%),  Vis,  Table  7)  denotes  the  change  in  classifier  accu¬ 
racy  between  the  first  quadrant  or  "quarter-coefficient" 
method  and  its  corresponding  "all-coefficient"  method.  For 
example,  60$  (+3%)  indicates  that  the  quarter-coefficient 
method  correctly  classified  3,%  more  samples  than  did  the 
corresponding  all-coefficient  classifier. 

When  using  (22),  to  classify  the  cloud  samples,  the 
quarter-coefficient  method  achieved  slightly  different  re¬ 
sults  than  did  the  all-coefficient  method.  There  were  no 
significant  changes  in  any  performance  rates  (±2  or  3  sam¬ 
ples),  suggesting  that  it  is  just  as  well  to  use  (22)  with 
power  spectra  generated  from  either  all  or  just  one-fourth 
the  coefficients  in  the  Fourier  transforms  of  each  imagery 
sample . 

When  using  (20)  and  only  visible  data  to  classify  the 
cloud  samples,  the  quarter-coefficient  method  achieved  high¬ 
er  results  (+7  to  +12  samples)  than  did  the  all-coefficient 
method.  Combined  visible  and  infrared  classifier  performance 
also  increased,  but  at  a  more  modest  rate  (+4  samples). 

The  results  for  this  data  set  favor  the  use  of  (20)  with 
the  power  spectra  generated  from  just  one-fourth  the  coeffi- 


TABLE  7-  NINE-CATEGORY  SPECTRAL  CLASSIFIER  ACCURACY  WITHOUT 


A 

PRIORI  AND 

COVARIANCE 

TERMS 

(22) 

" QUARTER -COEFFICIENT" 
(Compare  with  Table 

TRIAL 

3) 

Cloud  Type 

Number  of 
Samples 

Visible 

IR 

Combined 
Vis  and  IR 

Sc 

7 

57# 

29# 

71# 

St 

10 

40# 

20# 

80# 

Cu 

24 

63# 

71# 

79# 

Clear 

19 

74# 

79# 

79# 

Ac 

9 

33# 

56# 

56# 

As/Ns 

5 

0 

0 

80# 

Ci 

36 

56# 

72# 

83# 

Cs 

19 

74# 

79# 

'SS. 

ON 

CO 

Cb 

14 

86# 

71# 

86# 

Total 

143 

60# (+3#) 

64#( +3#)  80#( -1#) 

TABLE  8.  NINE- 

-CATEGORY  SPECTRAL  CLASSIFIER  ACCURACY  WITH 

A  PRIORI  AND  COVARIANCE  TERMS  (20) 

"QUARTER-COEFFICIENT"  TRIAL 
(Compare  with  Table  4) 

Cloud  Type 

Number  of 
Samples 

Visible 

IR 

Combined 
Vis  and  IR 

Sc 

7 

43# 

43# 

71# 

St 

10 

70# 

70# 

100# 

Cu 

24 

75# 

46# 

79# 

Clear 

19 

95# 

95# 

ON 

00 

Ac 

9 

78# 

22# 

67# 

As/Ns 

5 

100# 

100# 

1 00# 

Ci 

36 

oo 

69# 

72# 

Cs 

19 

89# 

79# 

89# 

Cb 

14 

86# 

71# 

86# 

Total 

143 

7  6# ( +  8# ) 

6 7#(-2#)  82#(+3#) 

TABLE  9.  FIVE-CATEGORY  SPECTRAL  CLASSIFIER  ACCURACY  WITHOUT 
A  PRIORI  AND  COVARIANCE  TERMS  (22) 

"QUARTER -COEFFICIENT"  TRIAL 
(Compare  with  Table  5) 


Cloud  Group 

Number  of 
Samples 

V  isible 

IR 

Combined 
Vis  and  IR 

Low 

( Sc/St/Cu ) 

41 

66* 

59* 

83* 

Mid 

(Ac/As/Ns ) 

14 

21  %> 

36* 

64* 

High 

(Ci/Cs) 

55 

73* 

00 

0 

Vi. 

ON 

00 

Clear 

19 

74$ 

79* 

79* 

Cb 

14 

86* 

71* 

CO 

On 

W 

Total 

143 

6?%(  +  l%o) 

69* (+2*) 

83*( -1*) 

TABLE  10.  FIVE-CATEGORY  SPECTRAL  CLASSIFIER  ACCURACY  WITH 
A  PRIORI  AND  COVARIANCE  TERMS  (20) 
"QUARTER-COEFFICIENT"  TRIAL 


(Compare  with  Table  6) 


Cloud  Group 

Number  of 
Samples 

Visible 

IR 

Combined 
Vis  and  IR 

Low 

(Sc/St/Cu) 

41 

71* 

71* 

90* 

Mid 

(Ac/As/Ns) 

14 

86* 

50* 

79* 

High 

(Ci/Cs) 

55 

00 

r\) 

78* 

91* 

Clear 

19 

95* 

95* 

89* 

Cb 

14 

86* 

71* 

86* 

Total 

143 

81*(+5*) 

75*(-2*) 

89*(+3*) 
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cients  of  the  Fourier  transforms  of  each  imagery  sample. 
However,  the  use  of  some  or  all  available  Fourier  coeffi¬ 
cients  should  not  be  confidently  determined  on  the  basis  of 
these  results  alone  since  the  data  sample  is  not  large.  As 
previously  mentioned,  time  savings  are  small.  Moreover,  no 
more  than  7%  of  the  samples  were  reclassified  as  something 
else  when  one-fourth  of  the  coefficients  were  used. 

4.7  The  Spectral  Classifier  and  3DNEPH 

The  discrete  Fourier  transform  spectral  analysis  al¬ 
gorithm  was  developed  as  a  potential  alternative  cloud  classi¬ 
fication  technique  to  that  part  of  3DNEPH  which  determines 
cloud  types.  Table  11  contains  the  spectral  classifier  (20) 
performance  rates  (all-coefficient  method)  for  those  cloud 
categories  used  by  3DNEPH.  Only  Cb  clouds  are  classified  by 
3DNEPH  when  using  visible  satellite  data.  Four  cloud  groups, 
Cs/Cb,  As/Ns,  Ci,  and  Ac,  are  used  when  only  infrared  data 
are  available  to  3DNEPH.  (A  slash  (/)  separating  two  cloud 
types  within  a  group  indicates  that  those  clouds  are  indistin¬ 
guishable  to  3DNEPH.)  The  eight  cloud  types  that  the  DFT 
classifier  uses  are  also  used  by  3DNEPH  when  both  visible  and 
infrared  data  are  utilized,  with  the  clear  case  omitted. 

On  the  average,  the  DFT  spectral  classifier  (20)  was 
able  to  correctly  identify  7-8  cloud  samples  out  of  every 
10  it  was  allowed  to  process. 

5.  CONCLUSIONS 

The  automated  DFT  spectral  classifier  was  allowed  to 
classify  143  cloud  imagery  samples.  The  results  obtained 
from  such  a  small  data  set  left  some  questions  unanswered. 
Nevertheless,  they  are  encouraging. 

The  classifier's  performance  rate  was  similar  for  two 
cloud  sample  subsets,  suggesting  that  it  is  consistent  in  its 
imagery  classification  decisions  and  that  the  classifier 
yields  reproducible  results. 


55 


— 1 


TABLE  11.  CLASSIFIER  ACCURACY  FOR  3DNEPH  CLOUD  CATEGORIES 


Cloud  Type 

Number  of 
Samples 

Number  of  Correct 
Identifications 

VISIBLE 

%  Correct 

CP 

14 

infrared 

11 

79% 

Cs/Cb 

33 

26 

79% 

Ci 

36 

26 

72% 

As/Ns 

5 

5 

100# 

Ac 

9 

4 

44# 

To  tal 

83 

COMBINED 
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VISIBLE  AND  INFRARED 

73# 

Sc 

7 

5 

71# 

St 

10 

10 

100# 

Cu 

24 

17 

71# 

Ac 

9 

7 

78# 

As/Ns 

5 

5 

100# 

Ci 

36 

23 

64# 

Cs 

19 

17 

89#* 

Cb 

14 

12 

86# 

Total 


124 


96 


77% 


By  partitioning  individual  cloud  types  into  larger 
categories  of  related  types,  many  of  the  classification 
errors  were  shown  to  be  classifications  to  related  cloud 
types  (e.g.,  Sc  for  Cu). 

The  use  of  infrared  spectra  to  determine  a  sample's 
cloud  type  proved  helpful.  When  used  alone,  the  IR-only 
classifier  correctly  classified  65-70$  of  the  imagery  sam¬ 
ples  in  the  data  set.  This  is  an  important  result,  since 
the  use  of  only  IR  data  is  forced  over  regions  underneath 
either  a  nighttime  satellite  pass  or  a  low  sun  angle  daytime 
pass.  As  hoped,  the  use  of  infrared  data  in  conjunction  with 
visible  data  greatly  enhanced  the  reliability  of  the  classi¬ 
fier. 

The  following  are  offered  as  suggestions  for  potential 
improvements  to  the  automated  cloud  imagery  classification 
algorithm : 

1.  The  averaged  2/3  nm  DMSP  data  used  by  this  routine 
provided  a  reasonable  spatial  resolution.  The  unaveraged 
1/3  nm  data  undoubtedly  contains  enhanced  small-scale  cloud 
feature  distinctions  and  should  be  tested;  however,  its  use 
may  not  yield  significantly  better  results.  The  use  of  1/3 
nm  imagery  would  require  the  routine  to  take  in  75  x  75- 
pixel  arrays  (25  x  25  nm  earth  coverage).  This  represents  a 
jump  from  37*37  =  1,369  to  75*75  -  5 >625  input  values,  an 
increase  of  over  four  times  as  much  data  that  would  have  to 
be  processed  by  the  routine.  The  cloud  identification  accu¬ 
racy  of  the  classifier  could  not  increase  fourfold  with  the 
implementation  of  1/3  nm  data,  while  data  processing  time 
would,  since  such  an  increase  in  the  amount  of  input  would 
place  additional  burden  on  computer  time  and  storage  allo¬ 
cations  . 

2.  The  collection  of  much  larger  imagery  sample  sets, 
adequately  representative  of  all  cloud  types,  would  firmly 
tost  the  oJassi  f'i  or’ s  ability  to  distinguish  among  cloud 
types . 

3.  The  many  background  features  that  show  through  the 
clouds  in  imagery  samples,  especially  those  containing  high, 
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thin  cirrus  or  fair  weather  cumulus  clouds,  cause  samples 
with  the  same  cloud  type  to  look  different  to  the  classifier 
if  their  backgrounds  are  different.  One  way  to  minimize  this 
adverse  effect  would  be  to  replace  all  cloudless  pixels  with 
a  brightness  value  of  0  (totally  black  pixel).  Samples  would 
then  have  a  common  background,  the  difference  among  all  sam¬ 
ples  being  only  the  cloud  types  in  them.  This  would  assure 
that  only  clouds  will  differentiate  one  sample's  cloud  type 
from  another’s. 

4.  The  spectral  classifier  had  a  performance  rate  nearly 
twice  that  of  the  mean  classifier,  but  the  spectral  classi¬ 
fier's  discriminant  function  required  26  times  more  informa¬ 
tion.  It  is  most  likely  that  not  all  26  components  of  the 
power  spectrum  of  a  3?  x  37  pixel  cloud  imagery  sample  need 

be  used  by  a  discriminant  function  to  obtain  similar  accuracy 
results.  Perhaps  adequate  discriminating  power  among  the 
cloud  types  lies  in  using  fewer  spectral  components,  the  only 
question  being  which  components  to  use  and  which  not  to  use. 

5.  The  main  feature  used  by  this  classification  al¬ 
gorithm  is  a  power  spectrum.  An  imagery  sample's  power  spec¬ 
trum  is  assumed  to  contain  enough  information  about  any  clouds 
within  the  sample  to  identify  the  sample's  cloud  type.  This 
was  obviously  not  always  true. 

It  was  observed  that  beyond  the  sixth  or  seventh  band 
number  p,  the  characteristic  spectra  of  all  cloud  types  had 
dissimilar,  distinct  ranges  of  values  (NAA^ax  -  NAA^n). 

When  compared  with  these  ranges,  a  sample's  corresponding 
observed  spectral  range  can  provide  additional  useful  infor¬ 
mation  to  the  automated  decision  process,  information  which 
can  successfully  identify  a  sample's  cloud  type  in  cases 
where  that  sample's  spectrum  is  not  sufficient  to  correctly 
classify  its  cloud  type. 

6.  As  wavenumbers  increase,  NAA’s  tend  to  decrease  in 
magnitude.  Consequently,  the  variances  of  the  higher  annu¬ 
lar  band  number  p  NAA's  also  tend  to  be  small.  In  this 
study,  generally  over  90%  of  the  variances  of  the  higher 
bandnumber  components  of  power  spectra  were  significantly 
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less  than  1 . 

Therefore,  the  productj^jjJ  (determinant  of  ^*1^)  of  all 
these  variances  (see  equation  (20))  was  an  extremely  small 
positive  number,  in  turn  making  the  term  -j  In  in  (20) 

rather  large.  It  was  often  greater  than  any  other  term  in 
(20),  including  In  P(w^).  Such  dominance  does  not  seem  to 
give  In  P(u»i)  the  weight  it  should  have  in  the  decision¬ 
making  process  of  the  classifier.  For  this  reason, 

"I  ln  PJ  is  better  left  out  of  (20)  or  scaled  down  to  a 
more  reasonable  magnitude.  The  a  priori  probabilities 
P(w  i) ,  since  £[0,l],  kept  the  term  ln  P(w^)  reasonable  in 
magnitude  and  should  remain  untouched  as  a  part  of  the  dis¬ 
criminant  function. 

7*  The  use  of  independent  data,  i.e.,  data  which  have 
not  been  used  to  calculate  and  (see  (16)),  would 

determine  how  sound  the  logic  of  the  classifier  is.  Mean 
spectra  and  covariance  matrices  can  be  generated  from  an 
imagery  sample  set  well  representative  of  all  cloud  types  in 
frequency  of  occurrence  and  in  texture.  These  spectra  could 
then  be  used  with  (20)  to  objectively  classify  different  sets 
of  imagery  samples,  and  would  eliminate  the  need  for  their 
subjective  classification  (unless  desired  for  independent 
verification  purposes).  In  its  current  state  the  classifier 
computes  |i^  and  directly  from  the  input  data  set,  in 

addition  to  requiring  interactively  supplied  input  in  the 
form  of  subjective  imagery  sample  cloud  type  identification. 

8.  In  this  study,  the  visible  and  infrared  spectra  of 
a  cloud  sample  were  used  simultaneously,  but  independently, 
to  identify  the  sample's  cloud  type.  The  use  of  cross¬ 
correlation  and  covariance  information  between  the  visible 
and  infrared  spectra  of  an  imagery  sample,  as  first  suggested 
by  Blackman  and  Pickett  (1979),  may  prove  to  be  a  sounder 
method  of  using  both  data  types.  Refer  to  Appendix  E  for  a 
description  of  this  method. 

On  the  basis  of  this  study's  results,  the  spectral 
classifier  has  established  itself  as  a  fast,  effective,  and 
consistent  satellite  imagery  classification  technique.  But 
as  with  any  prediction  technique,  this  one  will  ultimately 
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succeed  or  fail  depending  on  whether  or  not  the  results  it 
obtains  are  found  through  future  experiments  to  be  both 
accurate  and  reliable. 


i 
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Appendix  A 


RELATIONSHIP  BETWEEN  THE  COMPLEX  AND  TRIGONOMETRIC  FORMS  OF 

THE  DISCRETE  FOURIER  SERIES 


/A 

Let  F(i,j)  be  a  real-valued,  discrete  function  of  two 

variables  defined  on  an  M  x  N  region  of  M*N  equally  spaced 

grid  points,  where  M  and  N  are  positive  integers.  It  can 

/\ 

be  shown  that  the  complex  Fourier  series  of  F 


L2 


a 


=  E  E  *m.n' 


m=  [-"]+l  n=  [-f] 


r  N'  ,*mincc  . 


+1 


can  be  equivalently  expressed  in  trigonometric  form 
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=  2  f  +  2$R 

Tn.n  rm,-n 

=  -2  5^  +  2$  ^ 

Tn,n  Ym  , -n 

=  -2  -  2i|i  ^ 

Tn,n  *m,-n 

-  o  <R  .  o'T'  ^ 

’Hn.n  ^m,-n 


1  4  m  4  = 


“»4 


-  1 ,  and 


PROOF  OF  THE  PERIODICITY  OF  THE  DISCRETE  FOURIER  TRANSFORM 


It  is  not  immediately  apparent  from  equation  (6)  that 
the  discrete  Fourier  transform  of  the  periodic,  discrete, 
finite  function  F(i,  j)  (defined  in  Section  3.1)  of  two 
variables  is  itself  periodic,  discrete,  and  finite.  There 
are  only  M*N  distinct  complex  numbers  n  which  can  be  com¬ 


puted  by  (6).  Let  m  *  k,  n=  lbe  any  arbitrary  integers. 
From  (6):  M  N  ^  li 

$k ,  1  =  MN  A*  Z^F(l,j)e  M  N  . 


Now,  let  m  =  k+M,  n  =  1+N: 


*k+M,l+N  ~  M*N  4rff(i,j)e  L  M  N  J 


M  N 


m7n  Z  D(i,j)e’2,lA(kM  +  ^N} 


i=l  3=1 


e-2^i(Mi+Ni, 

M  N  .  . 

=  jji  E  E?(i,j)e'2,rA(1l  +  -N}e-2^1(i  +  J>. 
i=l  j=l 

Using  Euler's  formula,  and  since  i  +  j  is  always  an  integer, 
M  N  ... 

-  MR  +  “N5. 


1 


V  H 


v'" 


(cos  j)-  •£! 


i  +  j)) 


M  N 


Ef(i- j)e'2tr/rl(~l  +  ~i} 


i=l  j=l 


=  *k,l 


Therefore  there  are  only  M-N  distinct  numbers  for  which  (6) 

A 

can  be  evaluated.  Thus,  if  the  function  F(i,j)  defined  for 
an  M  x  N  grid  is  one  complete  cycle  of  a  discrete,  periodic 
function,  then  its  Fourier  transform  is  also  one  complete 
cycle  of  a  discrete,  periodic  function  with  M*N  unique  coef¬ 
ficients  in  each  period. 
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THE  DFT  OF  A  REAL-VALUED  FUNCTION 


In  order  for  the  complex  Fourier  series  expansion  of  a 
real-valued,  discrete,  finite  function  F ( i ,  j )  to  be  real¬ 
valued,  the  following  relationships  must  hold: 


di  =  I 

Tm,n  -m,-n 

\I>  -  ui 

rm,-n  r-m,n 

where  ♦  denotes  the  complex  conjugate  of  n-  These 

relationships  can  be  derived  by  expanding  (5)  and  setting  the 

coefficients  of  all  the  resultant  imaginary  terms  equal  to 

zero.  If  $ _  „  =  $  _  .  then  Pm  „  =  P  _  _  since  (from  (9)) 

m,n  -m, -n  ro,n  -m,-n 


•MJ*  <-£„>2 


li* R \2  T  /  *Im  71 
v  ^n,n  *  (Vn> 


=  P 

m,n 


65 


as&aattitt  m*  seats 


Thus  by  using  only  the  powers  of  those  coefficients  which 
lie  in  the  first  quadrant  of  the  m-n  plane  (as  does  the  "first 
quadrant"  method  described  in  the  Results  section),  any  infor¬ 
mation  contained  in  the  powers  of  those  coefficients  which  lie 
in  the  third  quadrant  is  used  as  well.  Essentially  one-half 
the  information  contained  in  the  transform  is  utilized  when 
one-fourth  the  available  coefficients  of  that  transform  are 
processed . 

Appendix  D 

CLASSIFIER  ACCURACY  MATRICES 

The  following  Figures  D-l  through  D-5  contain  matrices 
which  present  the  results  of  the  automated  classifiers  for 
the  set  of  143  cloud  imagery  samples  whose  composition  is  shown 
in  Figure  ?.  Each  figure  is  labeled  in  order  to  indicate  which 
classifier  variation  was  utilized  to  obtain  its  shown  results. 
The  cloud  sample  classifications  made  by  skilled  analysts  are 
represented  by  the  rows  of  each  matrix  (TRUTH),  and  the  deci¬ 
sions  made  by  the  automated  classifier  are  represented  by  the 
columns  of  each  matrix  (CLASSIFIER).  Agreement  between  the 
classifier  and  truth  is  shown  along  the  main  diagonal  (upper 
left  to  lower  right)  of  each  matrix.  Each  pair  of  matrices 
contains  a  frequency  distribution  on  the  left,  and  its  corres¬ 
ponding  relative  frequency  distribution  on  the  right.  The 
figures  show  in  detail  exactly  how  the  automated  classifiers 
performed . 

Overall  relative  frequency  for  all  or  any  group  of 
clouds  is  computed  by  summing  the  number  of  correct  hits 
within  the  cloud  type  or  cloud  group  of  interest,  and  then 
dividing  that  sum  by  the  total  number  of  samples  within  that 
cloud  type  or  group. 


VISIBLE  T 
R 
U 
T 
H 


CLASSIFIER 


Sc  St  Cu  u  Ac  <  Ci  Cs  Cb 


Sc 

2 

1 

2 

1 

1 

St 

1 

5 

2 

1 

1 

Cu 

3 

10 

9 

2 

CLEAR 

1 

6 

8 

2 

2 

Ac 

1 

3 

1 

1 

2 

1 

As/Ns 

1 

3 

1 

Ci 

8 

1 

6 

4 

3 

3 

9 

2 

Cs 

1 

1 

T 

1 

4 

9 

1 

Cb 

2 

3 

2 

3 

4 

CLASSIFIER 


Sc 

St 

Cu 

T  CLEAR 

R  Ac 
U  c 

T  As/Ns 
H  Ci 
Cs 
Cb 


a.  v> 

<  z 

Ui  \ 

Sc  St  Cu  o  Ac  <  Ci  Cs  Cb 


£ 

■ 

m 

■ 

■ 

£ 

m 

m 

ED 

m 

■ 

m 

■ 

■ 

m 

m 

m 

■ 

m 

m 

f*u 

■ 

■ 

m 

■i 

£ 

■ 

ed 

m 

ED 

ID 

m 

■ 

HD 

ESI 

m 

Q 

ID 

■ 

ID 

■ 

■ 

■ 

m 

■ 

£ 

m 

■ 

■ 

BS 

£ 

52 

SB 

ED 

52 

El 

m 

E 

■ 

123 

ID 

■ 

■ 

22 

m 

m 

52 

■ 

■ 

■ 

■ 

S3 

m 

m 

£ 

S3 

INFRARED 


CLASSIFIER 


Sc  St  Cu  d  Ac  <  Ci  Cs  Cb 


Sc 

2 

1 

1 

2 

1 

St 

2 

2 

1 

1 

4 

Cu 

2 

3 

14 

5 

CLEAR 

1 

4 

14 

Ac 

1 

2 

2 

2 

1 

1 

As/Ns 

1 

1 

2 

1 

Ci 

5 

5 

5 

4 

6 

6 

5 

Cs 

1 

1 

8 

9 

Cb 

1 

1 

1 

2 

9 

CLASSIFIER 

OC  at 
<  Z 
uj  \ 

Sc  St  Cu  o  Ac  <  Ci  Cs  Cb 
Sc 
St 
Cu 

T  CLEAR 

U  Ac 
T  As/Ns 
H 


Cs 

Cb 


E 

ID 

ID 

■ 

■ 

E 

ID 

■ 

■ 

E 

E 

■ 

IT*! 

■ 

m 

E 

■ 

■ 

52 

IE 

23 

m 

■ 

■ 

■ 

■ 

■ 

■ 

[2 

SJ 

ED 

■ 

■ 

■ 

■ 

■ 

ED 

£ 

■ 

m 

m 

o 

ED 

■ 

52 

■ 

■ 

ss 

E 

E 

■ 

■ 

E 

ED 

ID 

ED 

■ 

IQ 

IE 

■ 

ID 

■ 

■ 

52 

52 

■ 

■ 

■ 

EE 

m 

isa 

■ 

■ 

■ 

■ 

5D 

5S 

SD 

T 

CLASSIFIER 

CL  at 

<  z 

UJ  V. 

Sc  St  Cu  d  AC  <  Ci  Cs  Cb 


Sc 

2 

1 

1 

2 

f 

St 

1 

6 

1 

1 

1 

Cu 

2 

2 

14 

5 

1 

COMBINED  T  CLEAR 

7 

II 

1 

VIS  a  IR  R  A„ 

(J 

1 

2 

4 

1 

1 

T  As/Ns 

1 

1 

3 

Cl 

5 

3 

4 

3 

3 

15 

1 

2 

Cs 

2 

1 

1 

12 

3 

Cb 

1 

2 

2 

2 

7 

CLASSIFIER 


Sc 

St 

Cu 

T  CLEAR 

U  ^ 
T  As/Ns 
H 

Cl 

Cs 

Cb 


I  <0 

Sc  St  Cu  d  Ac  <  Ci  Cs  Cb 


IB? 

ID 

IH 

E 

ED 

■ 

■l 

m 

E 

EE 

13 

m 

■ 

■ 

■ 

m 

52 

S3 

01 

5S 

■ 

* 

■ 

■ 

■ 

■ 

53 

S3 

■ 

■ 

52 

■ 

■ 

m 

m 

■ 

■ 

EE 

■ 

ID 

ED 

■ 

■ 

E 

■ 

■ 

£ 

E 

■ 

■ 

■ 

EEID 

■ 

52E3EES3E; 

■ 

ED 

■ 

■ 

■ 

52 

52 

m 

SBj 

■ 

£ 

■ 

■ 

■ 

ID 

ID 

ED 

m 

Figure  D-l .  Mean  Classifier 


CLASSIFIER 

nr  <n 

<  ? 

Ld  >* 

Sc  St  Cu  o  Ac  «  Cl  Cs  Cb 


VISIBLE 


Sc 

St 

Cu 

CLEAR 

Ac 

As/Ns 

Cl 

Cs 

Cb 


4 

T 

2 

1 

4 

3 

2 

3 

16 

3 

2 

14 

3 

2 

1 

2 

2 

3 

1 

1 

3 

1 

4 

3 

4 

1 

15 

5 

4 

1 

T 

14 

1 

i 

1 

II 

CLASSIFIER 

C  (A 

<  z 

UJ 

Sc  St  Cu  o  Ac  <  Ci  Cs  Cb 


Sc 

St 

Cu 

T  CLEAR 
U  Ac 

T  As/Ns 
H 

Ci 

Cs 

Cb 


m 

■ 

SB 

■ 

■ 

0 

■ 

■ 

E 

■ 

■ 

■ 

■ 

S3 

■ 

m\ 

SB 

■ 

m 

■ 

■ 

■ 

m 

■ 

52 

■ 

■ 

m 

E 

■ 

■ 

SB 

in 

■ 

m 

■ 

m 

■ 

m 

■ 

■ 

ml 

■ 

■ 

■ 

■ 

■ 

E 

2) 

■ 

E 

33 

52 

■ 

ID 

'O 

■ 

SB 

ID 

■ 

■ 

■ 

52 

■ 

■ 

SI 

IE 

91 

■ 

52 

■ 

1 

52 

152 

m 

CLASSIFIER 


CLASSIFIER 


DC 

< 

UJ 


Sc  St  Cu  o  Ac  <  Ci  Cs  Cb 


R 

u 

T 

H 


Sc 

B 

... 

B 

D 

B 

St 

D 

B 

■ 

B 

B 

B 

Cu 

B 

■ 

U 

■ 

1 

□ 

□ 

CLEAR 

D 

B 

E 

■ 

■ 

B 

Ac 

i 

■ 

D 

B 

■ 

B 

2 

As/Ns 

■ 

B 

■ 

D 

O 

Ci 

~3l 

0 

■ 

m 

B 

2 

Cs 

1 

■ 

■ 

■ 

2 

03 

Cb 

c 

E 

□ 

2 

B 

D 

II 

<r 

< 

UJ 


<n 

Z 


Sc  St  Cu  o  Ac  <  Ci  Cs  Cb 


Sc 
St 
Cu 
T  CLEAR 
Ac 
As/Ns 
Ci 
Cs 
Cb 


m 

■ 

53 

E 

■ 

■ 

■ 

53 

■ 

Ii3 

■ 

m 

■ 

53 

ES 

■ 

23 

■ 

a 

■ 

52 

■ 

IQ 

IQ 

■ 

S3 

■ 

ID 

m 

■ 

■ 

■ 

m 

■ 

m 

■ 

■ 

ID 

■ 

■ 

m\ 

■ 

■ 

■ 

23 

■ 

i 

m 

E3 

■ 

52 

■ 

BE 

9 

m 

BE 

BE 

m 

■ 

■ 

■ 

s 

■ 

m 

a 

91 

■ 

■ 

■ 

E 

52 

m 

CLASSIFIER 

OC  in 

<  Z 
UJ 

Sc  St  Cud  Ac  <  Cl  Cs  Cb 


CLASSIFIER 


Sc 

St 

Cu 

COMBINED  T  CLEAR 


vis  a  IR 


R 

U 

T 

H 


Ac 

As/Ns 

Ci 

Cs 

Cb 


6 

1 

1 

7 

2 

1 

19 

1 

1 

1 

1 

2 

14 

1 

2 

1 

1 

T* 

1 

5 

2 

1 

1 

30 

2 

2 

17 

1 

_l_ 

12 

CC 

< 

UJ 


z 


Sc  St  Cu  d  Ac  <  Cl  Cs  Cb 


Sc 
St 
Cu 
CLEARl 
Ac 
As /Ns 
Ci 
Cs 
Cb 


m 

■ 

■ 

■ 

9 

9 

E 

9 

■ 

153 

S3 

■ 

■ 

9 

9 

E 

9 

9 

■ 

52 

E 

9 

52 

9 

BE 

52 

S 

m 

■ 

■ 

E 

9 

9 

22 

ED 

9 

IED 

■ 

E3 

9 

m 

9 

gj 

9 

■ 

■ 

■ 

■ 

9 

9 

m 

9 

9 

8 

BE 

52 

52 

9 

9 

9 

52 

ED 

m 

9 

52 

52 

Figure  D-2 .  Spectral  classifier  without  a  priori  and  co- 
variance  terms,  ’’all-coefficient"  method. 
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Figure  D-3.  Spectral  classifier  with  a  priori  and  covariance 
terms,  "all-coefficient"  method. 
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Figure  D-4.  Spectral  classifier  without  a  priori  and  covar 
iance  terms,  "quarter-coefficient"  method. 
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Figure  D-5-  Spectral  classifier  with  a  priori  and  covariance 
terms,  "quarter-coefficient"  method. 


CLASSIFICATION  USING  A  TRI-DIAGONAL  COVARIANCE  MATRIX 


The  combined  visible/infrared  classifier  used  in  this 
study  computed  the  26-component  power  spectrum  for  the  visible 
and  the  IR  of  each  37  x  37-pixel  imagery  sample  it  classified. 
Pickett  and  Blackman  (1979)  suggested  a  different  method  of 
combining  the  two  channels. 

Instead  of  a  26-component  visible  (x)  and  infrared  power 
spectrum  (y).  this  method  requires  that  a  2 ’26  =  52-component 
visible/infrared  spectrum  z  be  used  as  input  to  the  discrimi¬ 
nant  function  (16).  The  first  component  of  this  spectrum  is 
the  mean  visible  value  NAAq1S  of  an  imagery  sample,  followed 
by  that  37  x  37  sample's  25  visible  spectral  components 

V  i  q  . 

NAA^  =  xn,  n  =  1 ,2,3, . . . ,25-  The  2?th  component  of  z  is 

TR 

the  mean  IR  value  NAAQ  of  the  sample,  followed  by  25  infra- 

TR 

red  spectral  components  NAA~;  *  y  ,  m  =  1 ,2 , 3  ,  •  •  .  >25 . 

m  rn  > 

The  original  26  x  26  diagonal  covariance  matrix  /  , ^ 
would  increase  in  size  to  a  52  x  52  tri-diagonal  covariance 
matrix ^  composed  of  four  smaller  26  x  26  submatrices,  all 
of  which  are  diagonal.  The  main  diagonals  of  the  upper  left 
and  lower  right  submatrices  consist  of  the  variances  of  the 
visible  and  the  variances  of  the  IR  spectral  components  whose 
annular  band  numbers  are  the  same.  'The  main  diagonals  of 
the  lower  left  and  the  upper  right  submatrices  consist  of  the 
covariances  of  the  visible  and  infrared  spectral  components 
whose  annular  band  numbers  are  the  same. 

The  52 -component  column  vector  z  and  the  larger  tri¬ 
diagonal  covariance  matrix  are  then  plugged  into  equation 
(16),  with  r  =  52.  Thus  for  a  37  x  37  cloud  imagery  sample 
input  array, 


m  mm  * 


—  — 

“  *■ 

NAAq1s 

Z1 

NAA ^ 1 s 

Z2 

NAA^is 

• 

Z3 

• 

• 

NAA^S 

• 

• 

z26 

naaJr 

= 

z27 

naa™ 

z28 

NAA*R 

• 

z29 

• 

• 

NAA25  _ 

• 

• 

.  Z52. 

Evis 

^vis,  IR 

^pvis,  IR 

Ef 

so  that  (from  (16)) 

di(1)  =  4(  *  -  -  ^f1" 2*  -  iln  |E1 


+  InP ( to ■ )  i  i  =  1 ,2,3, ...  ,9* 

i  » 

The  column  vector  p^  is  the  mean  vector  of  all  vectors  i 
generated  from  samples  containing  cloud  type  i.  The  deci¬ 
sion  rule  would  choose  class  k  if  >  d^,  for  all  i  /  k. 
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GLOSSARY  OF  TERMS 


DMSP 

DFT 

R 

R2 

C 

z 

ZxZ 


F(i,  j) 


FFT 

* 

Re(z ) , Im( z ) , z 


x 

y 

Sc 

St 

Cu 

Ac 

As/Ns 

Ci 

Cs 

Cb 

Ruler's  formula 

n 


Defense  Meteorological  Satellite  Program 
discrete  Fourier  transform 
the  set  of  all  real  numbers 

the  set  of  all  ordered  pairs  (x,y)  such  that 
x  and  y  are  real  numbers 
the  set  of  all  complex  numbers 
the  set  of  all  integers 
the  set  of  all  points  (i,j)  such  that 
i  and  j  are  integers 
for  every . . . 

...is  an  element  of... 

the  greatest  integer  <  x,  x  positive; 

[-x]  =  -[x] 

a  discrete,  finite  approximation  to  the  con¬ 
tinuous  function  F(x,y) 
fast  Fourier  transform 
Let  z  =  a  +  bi  be  a  complex  number,  i 
«  *£*1 .  Then: 

Re(z)  =  a 
Im  ( z )  =  b 

*  .  . 
z  =  a  -  bi  (complex  conjugate  of  z) 

power  spectrum  of  a  visible  imagery  sample 

power  spectrum  of  an  infrared  imagery  sample 

stratocumulus 

stratus 

cumulus 

altocumulus 

altostratus/nimbostratus 
c irru  s 

cirrostratus 

cumulonimbus 

i  0 

e  =  cos  0  +  isin  0  ,  i  =  /T. 
determinant  of  the  square  matrix  ]£ 


n 


>1 


% 
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